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ABSTRACT 

Using SDSS Stripe82 data we have obtained deep radial surface brightness profiles of 7 face-on to 
intermediate inclined late- type spirals down to /ir' ~ 30 mag arcsec"^. We do not find any evidence 
for a sharp cut-off of the light distribution of the disks but a smooth continuation into the stellar halos 
of galaxies. Stellar halos start to affect the surface brightness profiles of the galaxies at ^ir' ~28 mag 
arcsec"^, and at a radial distance of > 4 — 10 inner scale-lengths. We find that the light contribution 
from the stellar halo could be responsible of previous classification of surface brightness profiles as 
Type HI in late-type galaxies. In order to estimate the contribution of the stellar halo light to the 
total galaxy light, we carried out a Bulge/Disk/Stcllar Halo decomposition by simoultaneously fitting 
all components. The light contribution of the halo to the total galaxy light varies from ~ 1% to ~ 
5%, but in case of ongoing mergers, the halo light fraction can be as high as '-^ 10%, independently 
of the luminosities of the galaxies. We have also explored the integrated {g' — r') color of the stellar 
halo of our galaxies. We find {g' — r') colors ranging from ~ 0.4 to ^ 1.2. By confronting these colors 
with model predictions, we encounter problems to fit our very red colors onto stellar population grids 
with conventional IMFs. Very red halo colors can be attributed to stellar populations dominated by 
very low mass stars of low to intermediate metallicity produced by bottom-heavy IMFs. 
Subject headings: Galaxies: Evolution, Galaxies: Formation, Galaxies: Spiral, Galaxies: Structure 
Galaxies: Photometry, Galaxies: Stellar Halos 



1. INTRODUCTION 

The faint outer disk and stellar halo of spiral galaxies 
hold a plethora of information related to their formation 
and evolution. This is because, due to their relatively 
long star formation and dynamical time-scales, the fos- 
sil records left behind by the formation processes survive 
longer than in the inner parts of the galaxies. Conse- 
quently, the characteristics of these faint subcomponents 
must be closely connected to the evolutionary paths the 
galaxies might have taken. Studying the structure and 
stellar populations of the faint outer disks and stellar ha- 
los makes possible to test the predictions of the current 
cosmological paradigm in detail. 

The disks of spiral galaxies have been studied in 
great detail. Structural studies of nearby spiral galax- 
ies often concern the inn er, bright stru c ture, the bulge , 
bar(s), rings, etc (e.g . Ide Jond Il996t iGrahamI 120011 : 
iMacArthur et al.ll2003D . During the past decade, how- 
ever, our knowledge of the structure and stellar pop- 
ulations of the disks at faint surface brightness levels 
has grown abundantly. These studies pushed the sur- 
face brightness limit to ^r' ~ 27 mag arcsec~^being able 
to explore the faint outskirts of disks. Thanks to these 
efforts, not only we got to know tha t a large amount of 
stars is present (IPohlen et all [20021: lErwin et al.) 120051: 
Pohlen fc Truiillol 120061 PT06 herafter: lErwin et al.l 



Based on their behavior galaxy surface brightness pro- 
files are classified into three major classes termed as 
Type I, Type II, and Type HI. Type I galaxies follow 
pure exponential profiles. The Type II morphology rep- 
resents all galaxies with a "downbending break" : a dou- 
ble exponential profile with a steeper outer exponential 
(jPohlenet al.|[2002l: [PTOllErwin et al.ll2008D . This class 
is the revisit ed and extend ed version of the Freeman 
Type II class (|Freemanlll970D now includ ing the so-called 
truncations discovered bv lvan der KruitJ (ll979'l. Type HI 
galaxies show an "upbending break" , a double exponen- 
tial pro file with a les s inclin e d outer exponential discov- 
ered bv lErwin et"aLl (|2005[) . IPT06I showed that only a 
minority of the surface brightness profiles of late-type 
disks follow a single exponential decline. The majority 
of late-type disks, ^ 90%, exhibits so called breaks in the 
faint outer regions. Apart from the fact that these fea- 
tures can be linked to different morphological features 
(such as bars, rings, spiral arms, etc) inside the disk, 
the properties of breaks follow tight correlations with 
the host galaxy properties, like Hubble-type, total lumi- 
nosity or stellar mass, and disk scale-length. Type II 
phenomenon is not limited to the local univers e, but 



has also been found in g alaxies at higher redshift (iPerez 
120041: iTruiillo fc Pohleni l2005). Moreover. lAzzollini et al 



200l " in these regions where star formation is negligi- 



ble, but they al so showed that t he paradigm of 'disk- 
as-expoiientiar (jPattersonl 119401 : Ide VaucouleursI 119581 : 
lFreemanlll970f) is far from accurate. 

|jbakos@iac.es, trujillo@iac.es| 



(I2008bir showed that the break radius could be used as 
a tracer of size evolution of disk galaxies. In the last 
~ 8 Gyrs, Type II disks went through significant (size) 
evolution. 

The stellar population properties obtained via 
color profiles revealed further clues and ca n be 
used to constrain models ()Roskar et al.l 120081 : 



Sanchez-Blazquez et alj l2009t iMartinez-Serrano et alj 
20091) that besides the effect of the motion of the 
stars (e.g. radial migration) consider the effects of 
star format i on, a nd/ or the thresh o ld in star formation 
(jKennicuttl Il989f ). iBakos et al.l ^00§) showed that 
these three surface brightness types directly translate 
into distinct behavior on their (g'-r') color profiles. In 
Type II galaxies the color profile follows a prominent 
"U-shape", with the bluest color at the break radius. 
This suggests that in Type II galaxies the origin of the 
brea k phen omena is related to star formation thresholds 
(e.g. ISchavd 120041 ). By having extended the analysis 
to intermediate redshifts, it is seen that the "U-shape" 
of th e (g'-r') color profiles is present in the last ^8 
Gyrs (jAzzoUini et al.ll2008aD . revealing that not just the 
morphology but the distribution of stellar population 
content is a common feature during the evolution of 
Type II galaxies. 

The surface brightness levels that are requiered to ex- 
plore the distribution of stars inside the disks are rela- 
tively bright {fir' ^ 27 mag arcsec"^), making it pos- 
sible to carry out photometric studies of disks to large 
distances (even to intermediate redshifts) on statistically 
representative samples. The stellar halos of spiral galax- 
ies are, however, so faint that to observe them it is nec- 
essary to go '-^ 10 magnitudes below the night sky level 
(to reach ~ 30 mag arcsec"^). This requires nothing less 
but excellent observing conditions free of observational 
artifacts such as flat-fielding problems or large skysub- 
traction uncertainties. Hitherto only a few integrated 
photometric studies can be found in the literature that 
have attempte d to reach the surface brightness of stel- 
lar halos (e.g. iZibetti et all 120041 : iJablonka et"all 120101 
and references therein). To investigate the statistical 
properties of s t ellar h alos [zlbctti & Ferguson (2004) and 
IBergvall et all ()2010D stacked ~ 1000 galaxies arriving at 
the level of fj,r' ^^-^31 mag arcsec"'^. These observations 
have also revealed anomalous red colors in stellar halos 
(otherwise known as the "red halo phenomenon") that 
could indicate the presence of extreme stellar populations 
(see also e.g. iLeaueux et al.l Il998t IZibetti &: Ferguson! 
[200l IZackrisson et al.ll2006l 120121) . 

Another way to explore the stellar population prop- 
erties of stellar halos is achieved by 'resolved star' 
technique. By prob ing, for instance, RGB stars inside 
the stellar halos (e.g . Ferguson et al.ll2 002: Davidge"2003'; 
Ferguson et al.l 12007: Reikuba et al.l i2009: Barker et al.. 
20121 and references therein), this method proved to be 
efficient in reaching surface brightness levels that are > 
2 magnitude deeper than what was so far accomplished 
by integrated photometry. However, the number of 
systems that can be studied by this technique is largely 
limited by the resolving power of current telescopes 
that only allows us to resolve stars in galaxies in our 
close vicinity (~ 15 Mpc). Despite of this, these studies 
provided stellar population and spatial contraints of 
stellar halos with unprecedented certainty. With the 
current cosmological paradigm in accordance, there is 
a growing consensus that the stellar halos of nearby 
galaxies such as M31, NGC 2403, NGC 0891 or M81 
are dominated by old but metal-poor populations. 
There is also compelling evidence that stellar halos 
are full of tidal debris, remnants of recent mergers 



with satellites of these giant galaxies. One of the most 
well studied galaxies, our neighbor, M31 went through 
quite a exciting merger history and hence is full of sub- 
structures (e.g. Ferguson ct al. 2002; Ib ata et al.l 120071 : 
IMcConnachie et al.ll2009l: iTanaka et al.ll2010D. similarly 
to th e stella r halo of t h e Milk y Way (e.g. [Ibata et al.l 
19941: iMa icwski c t al.l 119991: iNew bcrg ct al. 2001 



Gilmore et al.. .20021: 



Belokurov et al.l i2007i : .Juric et all 



200ilMIetaij 120081). There is also plenty of observa 



tional evidence showing ongoing merger phenomenology 
in spiral galaxies at these surface brightness levels, 
stellar streams have been detec ted around disk galaxies 
like NGC 5907, NGC 40 13 (IMartfnez-Delgado eFall 
I2008ll2009l) and NGC 0891 (|Mouhcine et al.ll2010D . 

One of the expectations of the current cosmolog- 
ical paradigm is that properties of the stellar halo 
(e.g. metallicity ) correlates with the properties of the 
host galaxv (e.g. |Brook et al.ll2004l : IBullock fc Johnston! 
120051: ICooper et al.ll2010[ ). Although, these correlations 
might be obscured due to stochastic variations in the 
merger/accretion history. Disks are expected to be less 
changed by external phenomena but more driven by sec- 
ular processes. To be able to distinguish between the 
consequences of these two entirely different effects, and to 
be able to discern the dominant mechanisms that shape 
spiral galaxies, we need to study the disks and stellar 
halos simoultaneously in as many systems as possible by 
reaching low surface brightness levels. 

The aim of this paper is to study the structure and 
stellar content of 7 late-type, intermediate inclined to 
face-on spiral galaxies down to a surface brightness level 
of fj,r' ~ 30 mag arcsec"^ using integrated light photom- 
etry that beforehand could only be reached by resolved 
star technique on individual galaxies. With this work we 
would like to demonstrate that with a correct data reduc- 
tion it is possible to study faint systems like the stellar 
halo around external galaxies as far as ~ 100 Mpc using 
the Stripe82 dataset. This opens the possibility to ex- 
plore a large variety of galaxy types and masses that are 
not present in our vicinity. 

2. THE SAMPLE 

We have created a sample of late-type spiral galax- 
ies selected from the HyperledcQ online catalog apply- 
ing several criteria for the Hubble-type (T parameter), 
the axis ratio (major axis/minor axis), the absolute B- 
magnitude (Mab^), and the radial velocity relative to the 
Local Group (corrected for virgocentric inflow, Vyir). 

We restricted the Hubble-type to cover the types be- 
tween Sbc and Sdm (2.99 < T < 8.49), building up 
an intermediate- to late- type sample. The axis ratio is 
selected to be logr25 < 0.301 (corresponding to an in- 
clination of < 61°). This allows us to have a face-on 
to intermediate inclined galaxy sample, and to mitigate 
problems related to dust extinction. The distance of our 
galaxies was restricted accordingly to create a volume- 
limited but co mplete sample of galaxies with Mats < - 
18.4 B-mag fsee lPTOel ). To this end, we selected galaxies 
that conform a radial velocity (corrected for the Virgo 
infall) Vvir < 3250 km sec~^ (out to ^ 46 Mpc esti- 
mated following the Hubble relation with Hq = 70 km s^^ 
Mpc~-^). Hyperleda yields a sample ~ 800 galaxies with 

http://lcda.univ-lyonl.fr/ 
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Table 1 

The Sample 



Gr&lclXy" 


(y. 





jnuDDie 


± 


i-J 


i\/r , 






(J2000.0) 


type 




[Mpc] 


[B-mag] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC0450 


01 15 30.4 


-00 51 40 


SABc 


5.8±0.5 


24.4 


-19.78 


NGC0941 


02 28 27.8 


-01 09 06 


SABc 


5.3±0.7 


21.9 


-19.13 


NGC1068 


02 42 40.7 


-00 00 48 


Sb 


3.0±0.3 


15.3 


-21.50 


NGC1087 


02 46 25.1 


-00 29 55 


SABc 


5.2±0.8 


20.7 


-20.65 


NGC7716 


23 36 31.4 


-1-00 17 50 


SBb 


3±0.5 


36.5 


-20.31 


UGC02081 


02 36 00.9 


-f 00 25 13 


SABc 


5.8±0.5 


36.5 


-18.53 


UGC02311 


02 49 27.9 


-00 52 23 


SBb 


3.0±0.4 


102.8 


-21.50 



no coordinate restrictions applied. To identify the galax- 
ies that fall within the area (~ 270 square-degrees) of 
Stripe82, we carried out a cross correlation between the 
full Hyperleda sample and the Stripe82 catalog. We also 
added a luminosity extended sample into this cross cor- 
relation, containing late-type spirals brighter than Mabs 
< -20.5 B-mag by allowing Vmr < 7000 km sec~^. 

In total in this paper we explored the following 7 spiral 
galaxies: NGC 0450, NGC 0941, NGC 1068, NGC 1087, 
NGC 7716, UGC 02081, and UGC 02311. The sample 
properties can be found in Table [TJ 

3. DATA AND PROFILE EXTRACTION 

3.1. The Data 

Our imaging data come from the Sloan Digital Sky 
Survey Stripe82. This stripe covers about ~ 270 square 
degrees of the sky in the equatorial plane and has been 
observed in all SDSS bands multiple, often more than 
80, t imes as part of the SD SS Supernovae Legacy Sur- 
vey (jAbazaiian et al.l 120091 ) . Even though this partic- 
ular survey was carried out to study variable objects, 
this implies that after stacking the data is ~ 2 magni- 
tudes deeper than the regular SDSS imaging. This makes 
Stripe82 presently the largest deep sky survey in the op- 
tical regime. 

3.2. Stripe82 Stacking 

The SDSS archive provides coadded imaging of 
Stripe82, but only of the first season of obse rvations (the 
SDSS Stripe82 coaddition is explained in lAnnis et alj 
1201 If ). We carried out different testing on these coadds. 
By having obtained surface brightness profiles of our 
galaxies from these images, it became evident that their 
coaddition process did not provide imaging with suffi- 
cient quality for our purposes. For this reason, we devel- 
oped our own pipeline to carry out (1) the skysubtraction 
on the downloaded chips; (2) the quality check of all the 
available runs; and (3) statistical rejection of low quality 
runs and the stacking. The pipeline is fully automatic 
and is mostly written in IDL, but makes use of a few 
external programs, as well. 

All the SDSS observations can be downloaded from the 
SDSS Data Archive Server (SDSS DAS). Before down- 
loading any image from DAS, we need to decide the 
size of the stack mosaic that include enough area to ex- 
tract the surface brightness profiles of our galaxies and 
estimate any residual sky affecting those. A field-of- 
view(FOV) of ~10 times the D25 proved to be sufficiently 
large to carry out the profile extraction. To download the 



images, we extract several information from CasJob^l- 
SDSS images are identified by their run, rerun, camcol, 
and field numbers. Casjobs returns these numbers cor- 
responding to the position of our target galaxies within 
a given search radius. (The search radius due to simple 
geometry is ~ V2xFOV/2) The downloaded Stripe82 
chips are flatfielded and bias-free, but they lack calibra- 
tion and still contain background. In our pipeline we did 
not include any calibration process, since what we need is 
a sufficient number of good quality imaging. We further 
explain the photometric calibration of the final stacks in 
Section 13.31 

After removing the 1000 counts of the SOFTBIAS 
added to all SDSS chips, the skyvalue is estimated in 
the following way. We measure the fluxes in about ten 
thousand randomly placed, 5 pixel wide apertures. We 
apply a resistant mean to the distribution of the aper- 
ture fluxes, and carry out several iterations to get rid 
of the fluxes biased by stars, and other background ob- 
jects. The mean of the bias-free distribution provides ac- 
curate measurement of the sky background, and is then 
substracted from the chips. The sky level and the back- 
ground scatter also serve as parameters to take into ac- 
count when assessing the quality of the imaging run. 

After obtaining the mean sky value and scatter of all 
runs, we study the distribution of these values. To avoid 
stacking runs with high background, for instance runs 
observed around full moon, we reject runs where the 
skylevel is a 3-sigma outlier from the robust mean of all 
runs. We also need to be careful about not to stack chips 
with relevant sky gradient, or chips which show signs of 
detector degradation. To assess skygradients and possi- 
ble chip anomalies, we use the same scatter information. 
In case of flat background the skylevel-skysigma relation 
is expected to conform the same relation for all the runs. 
For this reason, we fit this skylevel-skysigma relation and 
reject runs where it lies beyond 1-sigma error of the ex- 
pected scatter value. 

We also would have liked to avoid problems caused by 
wide-spread PSFs on our flnal stacks, to this end we only 
stacked runs with FWHM < 1.7 arcseconds. To this end, 
we measured the seeing on ever y run. We used a program 
called PSFEX (jBertinl [20T0al ) to extract the FWHM. 
This program makes a robust flt of th e PSF by using the 
stars extracted by SExtractor (Berti n fc Arno uts 1996). 
It is also important to avoid stacking of imaging with 
high extinction. Since Stripe82 observations do not come 
with photometric calibration, we needed to estimate the 

^ http:/ /casjobs. sdss.org/Cas Jobs/ 
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Table 2 





Stack and profile 


quafity 






Galaxy 


r'-band layers T^xp 


(min) 




l-^zLa — limit 


NGC0450 


(31)42 ~ 


38 


29.5 


28.8 


NGC0941 


47 ~ 


42 


29.9 


29.4 


NGC1068 


60 ~ 


54 


30.1 


29.5 


NGC1087 


55 ~ 


45 


29.0 


28.6 


NGC7716 


(38)58 ~ 


52 


29.4 


28.9 


UGC02081 


(43)53 


48 


29.8 


29.2 


UGC02311 


66 ~ 


60 


29.4 


28.9 



level of atmospheric extinction in a crude way. We ex- 
tracted stars from every run, and measured the fluxes in 
an aperture both on the Stripe82 run and on the DR7 im- 
age. The mean ratio of the measurements computed for 
all stars in the run was used as a first-order estimate of 
the extinction; we rejected a particular run, if we found 
more than 15% difference of the measured fluxes com- 
pared to DR7. Based their histogram, we found that the 
majority of the runs {> 60%) is within this flux differ- 
ence from DR7, hence this enabled us to create stacks 
using most of the observations. 

Since there have been Stripe82 observations taken 
under non-photometric conditions, this kind of quality 
check is necessary to have a final stack of exceptional 
quality. Once the non-useful runs were rejected, we start 
the stacking process which is done in two steps. During 
the first step, we create shallow mosaics of equal size to 
the final mosaic from all the fields contained within a 
given run (or more in case we need to use several cam- 
era columns to cover the desired field-of-view). The mo- 
saics are assambled from the skysubtracted chips using 
SWARP (|Bertinll20ToO ). Afterwards, we run SWARP 
over these shallow mosaics, applying a median stacking 
to all the available layers. It is important to note that 
the runs normally have slightly different offsets. Thanks 
to this, the scattered light from the telescope does not 
produce any pattern on the images either. Our final 
stacks do not show any residual sky-gradient or visible 
sky pattern. Due to the different number of layers used 
to compose the stacks, the mosaics have slightly different 
depths. Tabic [2] shows the number of layers used for the 
r'-band stacks (in brackets the number of g'-band layers, 
if that differs by more than 10%). For some of our galax- 
ies, the u'-band and z'-band detector problems prevented 
us to create reliable mosaics. In these cases (NGC 1068, 
NGC 1087, UGC 02081, and UGC 02311 in the u'-band 
and NGC 1068 in the z'-band ) we did not use the u'- 
and z'-band stacks for further analysis. 

3.3. Photometric calibration 

The photometric calibration of our deep mosaics is 
done by means of comparative photometry. The deep 
mosaics have an effective exposure time ~ equal to that 
of SDSS DR7. This makes it straightforward to compare 
number counts or surface brightness profiles to the DR7 
calibrated data. 

The SDSS DR7 photometric calibration is reliable at 
the level of '--^ 2 %. We opted to follow the traditional 
recipe to determine the zeropoint of our stacks by using 
the aa, kk, and airmass (the photometric zeropoint, the 
extinction coefficient and the airmass coefficients) . These 



parameters are provided in the tsField table associated 
with the chip containing the target object. The surface 
brightness zeropoint can be given as —2.5 x (0.4 x [aa + 
kk X airmass] + 2.5 x log{texp x pixelscale^), where the 
pixelscale is 0.396" /pixel and the exposure time for each 
pixel equals to 53.907456 seconds. 

We take the reliable part of the surface brightness pro- 
files of our galaxies obtained from the DR7 imaging and 
match it to the corresponding part of the Stripe82 sur- 
face brightness profiles. The mean difference is taken as a 
'photometric' offset and applied onto the DR7 zeropoint: 
ZPss2,x — ZPoR7,x + ^P-x- Since the the calibration 
is done in the same photometric filter system, no color 
term needs to be included. 

3.4. Masking and Profile Extraction 

We extract radial surface brightness profiles on masked 
mosaics in order to avoid contamination on our light pro- 
files. We apply conservative masking onto sources which 
clearly do not belong to the galaxy, like foreground stars, 
background galaxies, etc. These sources are extracted by 
SExtractor from a master image composed of all 5 bands. 
This master image is a stack of the different filters scaled 
to the r-band flux. Sources are extracted in "cold mode" , 
optimized for brighter objects. We use some SExtactor 
parameters such as the measured flux, elongation, and 
similar to determine the shape and size of these mask 
regions. Faint objects (fainter than ~ 20 mag) are not 
masked to statistically correct for the faint foreground 
stars unresolved from the galaxy image. In some cases, 
problematic sources, saturated stars, for instance, are 
masked manually. Since the masking is based on infor- 
mation coming from all 5 filters, we use the resulting 
mask as a mastermask, and it is applied to all filters. 

In order to extract radial surface brightness profiles 
corresponding to the underlying outer disk, we need to 
get representative values of the ellipticity and position 
angle of this region. Taking a s imilar approach to SEx- 
tractor (|Bertin fc Arnoutill996f) . we do this by comput- 
ing the second-order moment of the light distribution of 
the galaxy using the DR7 r'-band image. The second- 
order moment is directly related to the position angle, 
the semi-major (A) and semi-minor (B) axis lengths. 
We fix this ellipticity and position angle for all elliptical 
apertures. In order to increase the signal-to- noise coun- 
teracting its natural decrement due to the lower surface 
brightness in the outer regime, we used a logarithmic ra- 
dial sampling with steps of 0.03. In each aperture we 
estimate the galaxy fiux by the 3-(t rejected mean of the 
pixel values of that aperture. This helps to minimize the 
effect of morphological features like a spiral arm crossing 
the aperture. 

Our light profile, however, can still be contaminated 
by some residual sky background. This background is 
estimated by using equally spaced apertures. We obtain 
the number count profile of the galaxy up to very far 
distances, and we chose a large aperture where the pro- 
file becomes flat, beyond the identifiable profile of the 
galaxy. The S-cr rejected mean of the fiuxes inside this 
aperture gives a robust estimate of t he resid ual sky back- 
ground (for further information, see IPT061 their Section 
3.3). The error of the background determination is key 
to decide the surface brightness level down to which we 
trust our profiles. Following a conservative approach. 
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Figure 1. Left: RGB image of NGC 0941 created using the standard (single exposure) SDSS DR7 g'-,r'-,and i'-band imaging. Right: 
RGB image of NGC 0941 created using the g'-,r'-,and i'-band stacks of Stripe82 exposures of NGC 0941. The Stripe82 imaging is ~ 2 
magnitudes deeper than the standard SDSS imaging, this effectively means that S/N ~ 1 on the SDSS DR7 imaging corresponds to features 
of a surface brightness equal to 25 mag arcsec"^, however the same signal-to-noise on Stripe82 imaging corresponds to 27 mag arcsec"'^. 



this is placed where the profiles obtained by either over- 
or undersubtracting the sky measurement by icr start to 
deviate with more than 0.2 mag from the original profile. 
To see the limiting surface brightness of our profiles in 
the r'-band see Table [2] 

3.5. Data quality 

After stacking, Stripe82 imaging is expected to be ~ 
2 magnitudes deeper than the standard SDSS imaging. 
This amount of difference can be observed directly on 
the imaging: spiral galaxies show more outer structure, 
for instance, revealing outer disks extending farther out, 
meanwhile the inner morphology remains unchanged. As 
an example, see the DR7 and Stripe82 RGB color com- 
posits of NGC 0941 on Figured] 

A signal-to-noise ~ 1 per pixel on the DR7 imaging 
corresponds to features with a surface brightness of ~ 
25 mag arcsec"^. However, by summing up thousands 
of pixels in the outer apertures of the surface brightness 
profiles, we can reach S/N '-^ 5 at the level of ~ 27 mag 
arcsec"^. This is the S/N required to assure errors < 
0.2 mag. The same S/N ~ 1 per pixel on the Stripe82 
imaging can be associated with the surface brightness 
le vel of ^ 27 mag arcsec"^. 

IPT06I showed that with standard SDSS imaging it is 
possible to extract radial surface brightness profiles down 
to /ir' ^ 27 mag arcsec"^, hence we expect to be able to 
obtain reliable profiles down to /i^' ~ 29 mag arcsec"^. 
Table [5] shows that this depth is fairly stable for all our 
galaxies, however, in case of NGC 1068 or NGC 0941, due 
to a larger number of pixels falling in the outer apertures 
or smaller intrinsic noise thanks to better sky conditions, 
we were able to extract its profile down to /i^' ^ 30 mag 
arcsec"^. 

These deep surface brightness profiles also prove the 
high quality of the regular SDSS imaging. We made 
a direct comparison between the r'-band profiles of 
NGC 0941 in Figure H 



4. ANALYSIS 

4.1. Surface brightness profiles 

We obtained surface brightness profiles in all the five 
SDSS filters. These profiles are ~ 2 magnitudes deeper 
than regular SDSS profiles, giving way to explore the stel- 
lar populations of disks down to ~ 0.5 Mq pc~^. The 
deepest profiles, however, are obtained in the g'- and r'- 
filters, this is due to the better detector/filter response of 
these filters which results in smaller instrinsic noise. Due 
to this we explored in detail the r'-band surface bright- 
ness profiles and the {g' — r') color. Our results can be 
found in the Galaxy Atlas for every galaxy separately 
(see Appendix E]). 

4.1.1. Classification 

We used our deep r'-band surface brightness profiles to 
classify the surface brightness pro files of our galaxies. We 
follow ed the scheme presented in IPT06I and lErwin et al.l 
(|2008f) . The surface brightness profiles of Type I galax- 
ies are described by one single exponential, meanwhile 
Type II and Type III profiles both exhibit a well defined 
feature called the break. After an exponential decline. 
Type II profiles are followed by a downbending-. Type III 
profiles with an upbending-exponential after this inflec- 
tion point. There exist several sub-classes of Type II pro- 
files, related to the possible origin of the break. We dis- 
tinguish the so called Outer Lindblad Resonance breaks 
(Type ITORL), for these galaxies the break in the surface 
brightness profile lies close to the resonance zone caused 
by the bar, at ~ 2 times the bar radius. In case the break 
cannot be associated with such a well-defined phenom- 
ena, it is classified as Classical Truncation (Type II-CT). 
We show the results of our classification in Table |4l 

4.1.2. Excess of light at fi^' ~ 28 mag arcsec~^ 

The surface brightness profiles of our seven galaxies can 
be fitted into the aforementioned basic scheme. However, 
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Figure 2. Comparison of tiie surface brightness profiles of NGC 0941 obtained from standard SDSS DR7 imaging and Stripe82 stacks. In 
red we plot the surface brightness profile obtained from the standard SDSS imaging, the Stripe82 profile is overplotted in green. The shaded 
regions indicate the 1 and 3cr under- and over-substracted profiles in both cases. The surface brightness limits of our profiles are placed 
where the Icr under- and over-substracted profiles start to deviate by more than 0.2 magnitude from the original profile. The DR7 limiting 
surface brightness (26.9 mag arcsec"^) is indicated by the red dashed line, the Stripe82 limiting surface brightness (29.4 mag arcsec"^) is 
indicated by the green dashed line. On the joint panels we plotted the deviation of the DR7 profiles from the Stripe82 profile against the 
surface brightness and radius, this deviation is well below 0.2 magnitudes. This is a further proof of the reliability of the standard SDSS 
imaging. Stripe82 PSF: we also draw the PSF profiles of the g'- and r'-band imaging scaled to the r'-band central flux of NGC 0941. We 
observed no anomalous color gradient of the PSF. 



we note that most of our galaxies exhibit some extra Hght 
above the outer disk expectations starting at ~ 28 mag 
arcsec"^. This extra hght foUows a flattened profile that 
we think is not associated with any morphological feature 
belonging to the outer disk. In fact, in this paper, we 
work under the hypothesis that this extra light is diffuse 
light from the stellar halo. This extra light appears and 
shows the same structure in the other deep SDSS bands 
(g' and i'), as well. 

To further check our hypothesis, we investigated the 
origin of these features. First, we had to be sure that no 
artifact due to observational errors could possibly affect 
the quality of our profiles. There are several known prob- 
lems related to integrated photometry like flat-fielding, 
skysubtraction, contribution by saturated stars or ex- 
tended PSFs which could cause artificial excess (or de- 
fect) on our light profiles. 

Flatfielding — SDSS imaging due to their drift-scan tech- 
nique used during the observations show no evidence of 
flatfielding issues. 

Skysubtraction — On our surface brightness profiles we 
already took into account the possible effect of an error 
in the estimation of the backgound. We drew shaded 



regions corresponding to the la and 3cr under- and over- 
subtracted profiles (see Appendix |^ . The turnup point 
of the stellar halo feature (typically at ^r' ^ 28 mag 
arcsec"^) is at a surface brightness level > 15cr above 
the skylevel. The fact that we also see the same behav- 
ior for all the galaxies and in different bands, reinforces 
the idea that this is a real feature. 

Saturated stars — There is only one galaxy, NGC7716, 
where there is significant contribution from a saturated 
star that makes the interpretation of our surface bright- 
ness profile harder. This, however, could be minimized 
by extensive masking over the galaxy disk itself. We as- 
sessed the size of the region affected by scattered light 
by obtaining a pseudo surface brightness profile of this 
saturated star. We used similar approach to check the 
effect of saturated stars in case of other galaxies. 

Scattered light from the bulge — Given the depth of our 
images, we might expect the effect of the PSF on the 
bulge (the sharpest feature) of our galaxies to have a 
significant contribution at low surface brightness levels 
in form of an extended wing. This could cause the bulge 
light to scatter far enough to affect the surface brightness 



7 



profiles in the outer regions. To check the shape of the 
PSF at low surface brightness levels, we need to model 
accurately the PSF of the stacks. We used again PSFex 
to create the models of the PSF of our stacks. These 
PSF models are usually ~ 12 magnitude deep and might 
extend over < 50 pixels. By scaling these PSF models to 
the central surface brightness of the galaxy (see Figure 
[2]), we have seen that the extention and depth of the PSFs 
is negligible and could not affect the outer disk profiles. 
(We also considered the possibility of the bulge emerging 
above the disk, as a cause to observe the upturns on 
these profiles. Some arguments have been made in recent 
years that Type III profiles, that exhibit an excess of 
light over the disk in the outer regions, can be associated 
with the bulge emerging above the disk. The features we 
observe on the surface brightness profiles are a bit akin to 
Type III behavior, although in late-type spirals it would 
be somewhat unexpected to find bulges reaching ~ 10 
kpc or farther out.) 

After checking the above issues, we are confident 
that the excess of light observed at fir' ^ 28 mag 
arcsec"^ is real. We compared the surface brightness 
levels of our galaxies to several nearby galaxies like 
Andromeda (M31), M33, M81, or NGC2403. These 
galaxies are known to have well studied faint struc- 
tures which were identified by several inde pendent groups 
using 'resolved star' techniqu e (e.g. [ Ferguson et al " 
■2002, 2007; Ibata ct al. 2007; McConnachic ct al. 200 ' 
Barker et al.l 120091 120121) . The base of our comparison 



The central intensity of the inner and outer disks are 
Io,i and Io,2, and Q{Rf,r) is the step function: 
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is that our galaxies have quite similar stellar masses. 
iCourteau et al.l ()201lD showed that the regime of stellar 
halos actually starts at the level of /Xf ~ 27 mag arcse c"^ 
(using the color prediction^ of iVazdekis et al.l 120101 for 
old and metal-poor stellar populations, this magnitude 
translates into /i^' ~27.5 mag arcsec"^) corroborating 
our assumption of observing the signs of stellar halos on 
our profiles. 

4.1.3. Deriving disk and stellar halo parameters 

In addition to the classification of the surface bright- 
ness profiles, we carried out a detailed structural analysis 
on the observed r'-band profile. To derive reliable param- 
eters of the disk and stellar halo, a Bulge/Disk/Stellar 
Halo decomposition was done by simoultaneously fitting 
each component. With this we also avoid overestimat- 
ing the light contribution of the different components, 
especially that of the stellar halo to the total galaxy 
light. We fitted a composite model built up from func- 
tions corresponding to the known or assumed shapes of 
the different morphological components observed on the 
profile. It has been known for some time that disks 
are well fitted with exponential profiles (|Pattersonlll940l : 
lFreemanlll970[ ). However, the majority of disks does not 
follow a single exponential profile, but rather exhibits 
light deficiency or in some cases excess i n the outer disks 
(|van der KruitI 119791: lErwin et al.l [20051 [20081 [PT0 6) af- 
ter a well identifiable point of the profile, which is often 
called the break radius. Due to these reasons, we modeled 
the disk by a double-exponential to take into account the 
existence of these breaks: 



I{R) = /o.ie«/''ie(i?,,) +/o,2e^*/"Hl - e(i?,,)). (1) 

^ See also on |http://miles.iac.es| 



■br) 



1 a R< Rbr 
if i? > Rbr. 



(2) 



The step function assures the double-exponential pro- 
file with hi and /i2 scalelengths before and after R^r break 
radius. 

The bulges of spi ral disks are modeled wit h an r^/" 
profile (jAndredakis. Peleti er~fc Balcellsl 119951 ). The ID 
bersic r i/«-profiles (fSersid 119681) can be written as fol- 
lows: 



I{R) — leSXp ■ 



l/n 

^ I -1 




(3) 



where Re is the effective radius, and Ig = I {Re)- The 
quantity b„ is a function of the shape parameter n, and 
is chosen so that the effective radius encloses half of the 
total luminosity, a good approximation can be given as 
bn - 2n - 0.324 (e.g. [Graham 2001). 

Stellar halo profiles obtained by star-counts method 
are often fitted by Sersic-profiles, as well, (see e.g. 
ICourteau et al.ll20lH) . 

Altogether, the composite profile has 10 free parame- 
ters that needs to be fitted. In addition to this, in the 
flat, outer region where the stellar halo dominates the 
profiles, the observational errors are playing a significant 
role. Due to this, and after visual inspection, we decided 
to minimize the uncertainty in determining the proper- 
ties of the stellar halo by fixing its shape with a Sersic 
index n = 0.5, that is a Gaussian profile. iTrujillo et al.l 
(2002) showed that the mass density profile of the r^^"^ 
family can be extremely well approximated by an ana- 
lytical expression (see their Eq. 7) , which yields an exact 
case for n = 0.5. By choosing this constraint therefore we 
ensure a moderately declining mass distribution, on the 
other hand allowing a higher value for the Sersic-index 
would mean that the stellar halo has a prominent "over 
the disk component" in the inner parts of the galaxy (see 
their Fig. 2). 

The profile fitting is performed with the MPFIT rou- 
tine in the Interactive Data language (IDL), this rou- 
tine is based on the Levenberg-Marquardt algorithm 
(jMarkwardt 200^ to carry out non-linear least-squares 
minimization. Although this fitting method is extremely 
sensitive to initial conditions, with a proper set of input 
parameters the algorithm gives excellent results. 

The decomposition of the light profile yields robust es- 
timation of the scales of different galaxy component, like 
disk scalelengths and halo effective radius, that can be 
directly compared to theoretical expectations, and it also 
becomes possible to compare these halos to that of the 
Milky Way or M31. In case of the Local Group galaxies 
it becomes possible to integrate the observed starcount 
profile out to large radii, obtaining a robust and model- 
free extimation on the contribution of the halo light to 
total galaxy light. In our case, we only trace the inner 
part of the stellar halo, hence we cannot estimate the 
halo contribution by simple integration. However, one 
advantage of using Sersic models (as it is shown e.g. in 
[Tiruiillo. Graham, fc Caon.2001.) is that the global prop- 
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erties, such as the total hght (Ly) and the mean effective 
surface brightness ((A*)e), can be given in the form of ex- 
act analytical formula based on the aforementioned three 
observables (n, Re, and /e): 

Lt = kLieRl (4) 

where fc^ is the "structural parameter" that depends 
on the Sersic index: 

1, 27rn , , 

and 

^ -2.5^0,^. (6) 

The evaluation of the total galaxy light is as follows: 
we integrate the observed profile until the last measured 
point. Outside this region (to take into account the stel- 
lar halo contribution) we extrapolate the behavior of the 
fit of stellar halo to infinity. The stellar halo light frac- 
tions of the 7 galaxies compared to the host galaxy lu- 
minosity are shown in Figure [H 

4.2. Stellar Population Analysis 
4.2.1. Color Profiles 

In addition to the light distribution of these late-type 
spirals we obtained color profiles to explore the radial 
distribution of stellar populations in our galaxies. Nor- 
mally, color profiles cannot be studied down to the same 
surface brightness level as the surface brightness profiles, 
due to the fact that errors propagate quadratically and 
the uncertainty of the color is necessarily attributed to 
the uncertainty of the sky estimation in both filters. The 
errors coming from the sky estimation, however, can be 
alleviated by using large bins. We find that the most 
convenient strategy is to place the color bins around the 
well identified morphological features such as the bulge, 
disk etc. Nevertheless, it is important to avoid using a 
single bin to describe regions where the color gradient 
is steep, and it would mean mixing very different stel- 
lar populations. The bins we applied for the g'-r' color 
profiles (see Appendix |Xl are a compromise between the 
increase of signal-to-noise and underlying color gradient. 

In the regime of the stellar halo we observe no signifi- 
cant color gradient, therefore to give a robust result, we 
obtained a single (g'-r') color for our stellar halos us- 
ing one large bin. To derive these colors for all stellar 
halos in a consistent way, we placed the inner bound- 
ary of the color bin at the radius (R>5o) where the stel- 
lar halo light starts to dominate the total light, and the 
outer boundary at the radius where the surface bright- 
ness profile reaches the 3cr skynoise. This way the binsize 
depends on the stellar halo properties and instrinsic noise 
properties of the image. 

We also obtained color profiles for other adjacent filters 
(u'-g',r'-i',i'-z'), in total providing 4 color profiles that 
can be studied in detail. 

The Galactic extinction has been taken into using the 
Schlegel et al. (1998) values. We used the IDL code pro- 
vided by Schlegel et al|3 to get the values corresponding 

http:/ /astro. berkeley.edu/~marc/dust/indcx. html 



to the different galactic positions of our sample. 

4.2.2. Color errors 

Deriving colors for the faint stellar halo calls for a lot 
of care. At very faint surface brightness levels, the col- 
ors are affected by the shape of the PSF of the different 
bands. Anomalous PSF scattering can cause excess of 
light in one band, hence caus ing an artificial color. This 
effect was presented e.g. by Ide Jong (2008) who accu- 
rately proved that (1) PSF tails on stacked SDSS imag- 
ing can extend as far as '-^ 100 pixels reaching ~ 14 mag- 
nitudes below the central surface brightness, (2) color 
gradients and reddening could be artificially caused by 
the extended scattering especially in the SDSS i'-band. 

Particularly , in t he case of edge - on g alaxies 
(|Zibetti et all 120041: iZibetti fc Fergusonl 120041 ) . ac- 
cording to the Ide Jonel (|2008l ) results, the scattered light 
from the central disk contaminates the light of the stellar 
halo. In face-on projections, however, the artificial color 
gradient induced by the different wing extentions of the 
PSF plays a less relevant role. This is for two reasons: 
(1) the slope of the surface brightness profiles in face-on 
projections are much moderate than in the edge-on case 
being less affected by the PSF; (2) whereas the stellar 
halos are found after the disk inmediately in edge-on 
projections, in face-on projections the transitions are 
much moderate and the PSF size is not large enough to 
affect the stellar halo in their outer regions. 

Colors could also be affected by an inaccurate sky sub- 
traction. To account for the propagation of the uncer- 
tainty of the sky estimation into the color, we over- and 
undersubtracted the galaxy images with la uncertainty, 
and obtained total intensity, hence magnitudes in the 
different bins from those images. The mean magnitude 
difference (6) = |m+o- - ni-o-l/S was propagated quadrat- 
ically: 

and used as error on the color. 

4.2.3. Stellar surface mass density profiles 

To study the stellar mass distribution of our spirals, 
we computed the ste llar mass density (E^) profiles using 
the iBell et all ()2003[ ) prescriptions for the mass-to-light 
(M/L) ratio. 

In this work we assum e a Kroup a IMF (jKroupal 
1200 ID . which according to IBell et all (j2003D implies a 
deduction of 0.15 dex from the mass-to-light ratio using 
the following expression: 

logw{M/L)x = (ax + bx x color) - 0.15, (7) 

where for the SDSS (g'-r') color ax = -0.306 and bx = 
1.097 is applied to determine the r' band M/L. 

Once we know the (M/L) ratio, using the expression 
below it is straightforward to link the stellar mass density 
(S*) profile with the surface brightness profile at a given 
wavelength (/^a)- 

logio^. = loglo{M/L)x-OA{^lx-mabs,Q.x)+8.629, (8) 

where mabs,Q,\ is the absolute magnitude of the Sun at 
wavelength A, and is measured in Mqpc~^. 
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Figure 3. Left panel: Fraction of the stellar halo light vs. the total galaxy light. The stellar halo light is estimated by fitting with a n=0.5 
Sersic-mode the region of the stellar halo of the galaxies, meanwhile to obtain the total galaxy light we integrated the observed profile. The 
dark red points coincide with the galaxies where we i dentify the presenc e of tid al features (NGC7716) or m inor mergers (NGC1087). The 
stellar halo light fraction observed in the Milky Way l lCarol'lo et al.l|20i a). M31 (Courteau ct al. 2011?), M33 fMcConnachi e et al.| [ 2010l ) and 
NGC2403 (Barker ct al. 2012!) are also overplotted. The light contribution of the stellar halo to the total galaxy light varies from ~ 1% to 
~ 5%, but in case of ongoing mergers the stellar halo light fraction can be as high as ~ 10%. Right panel: Effective surface brightness of 
the stellar halo vs. the integrated (g'-r') color of the stellar halo. When computing the stellar halo color, the halo region is taken outside 
the radius where the halo light start to contribute with > 50%. We find colors ranging (g' — r') ~ 0.4 - 1.2. The extremely red halo colors 
(g' — r' >0.8) arc found in galaxies with no evident interaction. 



4.3. Stellar population models 

Comparing the observed stellar halo color to the pre- 
dictions of single stellar population (SSP) models we ex- 
pect to get a crude hint on the stellar population that 
might build up the observed stellar halos. Since we lack a 
wide multiwavelength data, we opted for constraining the 
observed (g'-r') colors. We used the latest MILES SSP 
models. These empirical models are based on real stellar 
spectra and give e xcellent cove rage both for metallicity 
and ages (Vazdeki s et al.|[2010[ ). Another important as- 
pect of using these models is that they provide predic- 
tions for different IMFs. 

We have decided to constrain the stellar populations 
of the stellar halos in a novel way. This was necessary 
since the age-metallicity degeneracy prevents us to gain a 
deeper understanding of these stellar populations, more- 
over we would like to see how the observational errors af- 
fect the measurements of age and metallicity. To this end 
we created synthetic color maps based on the Vazdekis 
models. 

As a first step, we took all SSP model spectra and 
convolved them with the SDSS g' ad r' filters in order to 
compute model colors for each age and metallicity. Then 
we composed a two-dimensional map of the model colors 
and placed our measurements over the age-metallicity 
grid. To visualize how the age-metallicity degeneracy 
appears in the colors of the stellar population, we as- 
signed the colors of the rainbow table in IDL to the values 
computed from the models. The rainbow colors change 
smoothly from violet (black) to red. Instead of comput- 
ing values from the observed colors, we decided to 
use the observed colors and its errors, directly on these 
maps. The range of colors within the observed color and 



its error bar can be associated with a range of ages and 
metallicities, creating by default a two-dimensional prob- 
lem. This is shown in Figured 

We have found that some of our observed stellar halo 
colors are so red, that they lay outside the grid of mod- 
els based on the conventional Kroupa IMF with slope a 
= 1.3. The observed halo colors are redder than the 
reddest colors of those models. Redder model colors 
can b e expected of m odels based on bottom-heavy IMFs 
f Vaz dekis et al.|[2010( ). Although, there is a wide range of 
IMF slopes between a = [0.3, 3.3] available, here, we sim- 
ply took the model that provided us with the reddest col- 
ors, therefore we used models with Unimodal IMF with 
a slope of a = 3.3. To comprehend how the IMF affects 
the intrinsic color distribution, see Figure SI 

5. RESULTS 

5.1. Surface brightness profiles 

The surface brightess profiles of our 7 galaxies have 
been classified accordingly to the scheme explained in 
Section [4X11 

Our sample contains one Type I disk 
NGC 1068(=M77). In this work, NGC 1068 has 
been classified as Type I (in 1PTO6 it was classified as 
Type II-ORL). We considered here that rather than 
having a break as sociate d with the Outer Lindblad 
Resonance cited by IPT061 the extended oval distortion 
in the disk accounts for the extra light over the disk 
observed on the surface brightness profile (see Figure 
[TU]). NGC 0450, NGC 0941, and UGC 02081 exhibit 
Type II classical truncations, the break of UGC 02311 
lies by the resonance zone caused by a bar, hence it is 
classified as Type II-ORL. When classifying NGC 1087 
and NGC 7716, we had to be a bit cautious. These 
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Figure 4. Metallicity-Age maps; Using the MILES SSP models (Vazdekis et al."2010') to get a liint on what stellar populations might 
build up the stellar halos. We illustrate the expected age and metallicity for two different galaxies. In case of NGC 0941 the stellar halo 
is extremely red and models produced by regular Kroupa IMF are unable to reproduce this color. The blue (NGC 0941) and magenta 
(NGC 1068) dashed lines show the age and metallicity solutions that are compatible with the observed colors within their error bars. The 
upper panel displays the model colors for Kroupa IMF and the bottom panel with a bottom- heavy IMF. 
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Table 3 

Bulge Properties 



Galaxy 


n 




Me,!) 






[kpc] 


[mag arcsec"^] 


NGC0450 


0.95 


0.94 


20.77 


NGC0941 


0.64 


1.34 


21.73 


NGC1068 


0.92 


0.54 


16.54 


NGC1087 


0.73 


0.24 


17.81 


NGC7716 


1.50 


0.60 


18.18 


UGC02081 


0.50 


0.45 


20.49 


UGC02311 


1.33 


4.18 


20.21 



galaxies seem to have a break feature in the inner 
disk that can be associated with a centrally more 
concentrated star forming region. This however is not 
the kind of phenomena we designate as Type II, hence 
these galaxies are classified as Type II*. In their outer 
disk both galaxies exhibit Type III up-bending. In 
NGC 7716 this extended Type III feature has a disky 
origin starting at about 80 arcseconds at a level of 25 
mag arcsec"^, meanwhile in NGC 1087 the up-bending 
shows a spheroidal origin. To be consistent with pre- 
vious classification scheme, these 2 galaxies have been 
classified as Type II-I-III. 

Interesting to note that the high contribution of the 
stellar hal o light in case of NGC 1087 explains why this 
galaxy in IPT06I was classified as Type III. However, in 
this sample, we find that the presence of the stellar halo 
causes a "Type III like" feature on all profiles. 

Beyond the classification, we derived several properties 
related to the radial light distribution of these galaxies. 
Applying the decomposition explained in Section 14.1.31 
we inferred several structural properties for the bulge, 
disk, break, and stellar halo. 

5.1.1. Bulge properties 

This paper was not aimed to discuss bulges in de- 
tail, however, it is interesting to note that the bulges 
of our galaxies are well fitted with low Sersic-index(n) 
values n < 1.5 (see Table [3]). IKormendv fc Kennicutt] 
(|2004l ) discussed that small n values are signatures of 
pseudo-bulges. There is also a genera l trend noted by 
lAndredakis Peletier. fc Balcellsl (|1995f ) that in late-type 
(Sc and later) galaxies the bulge profiles are very close 
to being pure exponentials. We observe the same trend 
for our objects, and that indeed in some cases, like 
NGC 0941, we have alternatively confirmed the presence 
of a pseudo-bulge (Kormendy, private communication). 

5.1.2. Disk and Break properties 

Our model fitting provides robust estimations of disk 
properties, such as inner, and outer scale-lengths, break 
radius, inner and outer disk central surface brightness, 
and the break surface brightness. The results of the 
fitting can be found in Table SI A detailed statistical 
analysis of these parameters is beyond the scope of this 
paper, due to the low number of observed galaxies in our 
sample. However, on a one-by-one basis it is interesting 
to compare their properties to previous studies of late- 
type disks, and later on, to put this in context with the 
o bserved stellar halo properties. 

IPT06I studied ^ 90 galaxies where they have found 
that ~60% of late-type disks contain Type II breaks. 
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Figure 5. D isk property correlations with absolute magnitude 
for the IPT06I Classical Truncation sample. The red dots are the 
properties obtained for the Type II galaxies of the current sample. 
Overplotted in both panel are robust linear fits to guide the eye and 
in the upper right corner the Spearman rank correlation coefficient. 
Upper panel: The inner and outer disk scale-length ratio for the 
IPT 06 sample of Type II-CT galaxies. Bottom panel: The break 
radius in inner scale-length units. 

Based upon this result it is not surprising that four of our 
seven galaxies (-60%) exhibit Type II behavior. IPT06| 
also studied the typical radius where surface brightness 
breaks usually occur, taking into account that the mor- 
phologically distinct sub-classes might have very differ- 
ent origins. Five of our galaxies (NGC 0450, NGC 0941, 
NGC 1068, NGC 1087 and UGC 02081) were part of the 
PT06 sample, as well. Besides of the re-classification of 
NGC 1068 and NGC 1087, the break radn we measured 
for NGC 0941(R6^ - 72") and UGC 02081(Rbr ~ 67") 
are different from the PT06 values (Rtr = 82" and 53", 
respectively) . Part of these differences are due to the dif- 
ferent ellipticities we applied when obtaining the surface 
brightness profiles. Different ellipticity means different 
scale-length properties, as well. However, differences in 
the measurement of the scale-length might also arise due 
to a more accurate Bulge/Disk decomposition method 
we used to infer disk properties. We nevertheless are 
in position to say that statistically the obtained break 
properties are within the PT06 expected ranges. In light 
of both results, we conclude that the method applied in 
PT06 was sufficient to give contraints on the disk struc- 
ture on a statistically significant level, and that the use 
of shallow DR7 data does not compromise disk measure- 
ments. 

Similarly to lPTOl we characterized the strength of the 
break by the ratio of the inner and outer scale-lengths. 
The mean ratio of the inner and outer scale-length on 
the PT06 sample was 2.5±0.5. They also reported a 
lower range of values that could be as low as 1.1 (case of 
NGC7716). To compare how our measurements fit into 
the lPT06l relations, we overplotted the of this paper 
w ith the PT06 values (see the upper panel of Figure [5]). 

IPT06I showed that Type II Classical Truncations 
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Table 4 

Disk and Break Properties 



Galaxy 


SB-Type 
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NGC0450 


Type II-CT 


82.3 


9.7 


2.20 


37.4 


4.4 


2.21 


21.19 


18.29 


23.53 


0.26 


8.2 


NGC0941 


Type II-CT 


71.5 


7.6 


3.45 


20.7 


2.2 


1.72 


20.82 


18.13 


24.50 


0.35 


4.2 


NGC1068 


Type I 








71.5 


5.3 




20.91 










NGC1087 


Type II*+III 


51.4 


5.2 


1.37 


37.4 


3.8 


1.65 


19.77 


18.81 


21.17 


0.39 


98.8 


NGC7716 


Type II*-|-III 


39.2 


6.9 


2.53 


15.5 


2.7 


1.10 


19.85 


19.59 


22.47 


0.40 


30.8 


UGC02081 


Type II-CT 


66.5 


11.8 


3.36 


19.8 


3.5 


2.22 


21.60 


17.16 


25.32 


0.28 


1.7 


UGC02311 


Type II-ORL 


42.1 


21.0 


2.38 


17.7 


8.8 


2.30 


21.47 


18.11 


23.92 


0.43 


8.8 



(NGC 0450, NGC 0941 and UGC 02081) are expected 
to occur at '-^ 9 ± 2.4 kpc, meanwhile breaks related to 
resonances (UGC 02311), show a wider span of galac- 
tocentric distances with the same mean distance (~9 
kpc). Both distributions showed the tendency of larger 
break radius occuring in higher luminosity galaxies (see 
their Figure 7). Taking into account their results, it 
is not surprising to measure such a large break radius 
in UGC 02311, since this galaxy is ~2-3 magnitudes 
brighter than the other galaxies with Type II-CT breaks 
(see Table [T|). The break radius measured in relation 
to the scale-lengths, however, does not reveal such differ- 
ences between the two Type II sub- types, its typical value 
being 2-3 him this mach es well with the mean parameter 
value of 2.5±0.6 found in lPTOGi (see lower panel of Figure 
[5]). The inner scale-lengths of these galaxies co rrespond 
to the typical values measure d in local galaxies (jde Jonal 
\199^: 'Mac Art hur et al.l[200l . although UGC 02311 has 
a scale- length larger than usual (~ 9 kpc). 

The inner disk central surface brightnesses are ^ 0.5 
magnitude fainter compared to the PT06 values. This 
difference could be attributed to the better treatment of 
the light contribution from the bulge done in this pa- 
per. The wider spread deviations in the outer disk cen- 
tral surface brightnesses compared to the PT06 values 
are more difficult to explain. The proper fitting of the 
stellar halo should in any case cause the outer disk scale- 
length to shorthen and its central surface brightness to 
brighten, resulting in a sharper break. However, PTQ^ 
breaks appear to be sharper (having brighter central sur- 
face brightnesses, and larger hin/hout), probably the ef- 
fect of the ellipticity, that we intend to investigate in our 
upcoming paper (BT2012, in prep.). 

In order to investigate how Type II break proper- 
ties (namely the inner and outer scale-length ratio) de- 
pend on the observing wavelength, we carried out a 
Bulge/Disk decomposition on the profiles in all 5 filters 
(if available) by applying the bulge and di sk part of the 
fitting routine explained in Section 14.1.31 To be fully 
consistent, we fitted the same disk region for all filters 
for our galaxies. We have cut the profiles taking care of 
not mixing the regions where the stellar halo light would 
have a significant contribution. During the fit we fixed 
the break radius. 

We find that the hin/hout decreases smootly from u'- 
to z'-bands in our Type II galaxies (see Figure [6]). 



5.1.3. Stellar Halos 



Using SDSS Stripe82 data gives us the opportunity to 
study the stellar content of disk galaxies down to a very 
deep level ~30 mag arcsec"^ in the r'-band. This depth 
enables us to get a glimpse of the stellar halo on the 
surface brightness profiles: the light contribution from 
the stellar halo leaves a well identifiable upturn on the 
r'-band profiles, and when the stellar halo starts to dom- 
inate the light the profile becomes fiatter. We observe 
this flattening on the profiles of all 7 galaxies that helped 
us to establish that the typical surface brightness level at 
which stellar halos start to affect the shape of the surface 
brightness profiles is ^27 mag arcsec"^. 

Although, this is a purely observational concept, it is 
interesting to discuss the radius where the light profile of 
the (outer) disk of the galaxy first starts to deviate from 
the exponential decline. We refer to this radius as the 
turn-off radius (Rtumoff) in Tabled The estimation of 
this radius (and surface brightness) is model-dependent. 
Based on the assumption that the disk continues expo- 
nentially without having an actual edge, a 0.5 magnitude 
change (larger than the error bars of our measurements 
at these levels) could be expected if the halo light con- 
tribution reaches ~ 25% that of the disk. This 0.5 mag 
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Figure 6. The strength of the break (indicated with the inner- 
outer disk scale-length ratio) obtained from the surface brightness 
profiles observed in the 5 SDSS bands and from the stellar surface 
mass density profile (See l4.2r3ll . The inner-outer scale-length ratio 
decreases rapidly towards the redder filters, moreover on the stellar 
surface mass density profile the break disappears {hin/hout ~1). 
However, in case of UGC02081 the break remains prominent on 
the stellar surface mass density profile, implying that the break in 
this galaxy is entirely different phenomena and reseambles more of 
the cut-offs found by Ivan der Krui"3 l|197SI ). 
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Table 5 

Properties of stellar halos 



Galaxy 



[kpc] 



[mag arcsec"^] 



^^tumoff 



[kpc] 



[kpc] 



NGC0450 
NGC0941 
NGC1068 
NGC1087 
NGC7716 
UGC02081 
UGC02311 



18.90 
9.90 
38.00 
13.50 
35.00 
20.11 
40.00 



27.64 
28.14 
28.14 
25.46 
27.74 
28.72 
27.97 



128.68 
100.05 
393.82 
111.46 
86.96 
84.08 
62.98 



15.22 
10.62 
29.21 
11.19 
15.39 
14.88 
31.39 



3.44 
4.83 
5.51 
2.98 
5.61 
4.25 
3.56 



152.48 
121.20 
513.78 
159.11 
106.89 
97.36 
74.62 



18.04 
12.87 
38.11 
15.97 
18.91 
17.23 
37.19 



deviation from the expected exponential shape of the sur- 
face brightness profile is well observable. The measured 
values of Rtumoff fall at 10 to 30 kpcs from the center 
of the galaxy, in terms of the inner scale-length between 
3 and 5.5 hin. 

These upturns should not be confused with the outer 
regions where the stellar halo clearly dominates the light 
of the galaxy (since around the upturn the disk is still 
the major contributor of the galaxy light). The radius 
where more than 50% of the galaxy light comes from the 
stellar halo is referred as R>5o. At this radius the stellar 
halo contribution clearly overcomes that of the disk and 
dominates the flat appearance of the surface brightness 
profile. As it is seen in Table \5\ there are, on average, 
2-3 kpcs differences between Htumoff and R>5o- 

In all these aspects, NGC 1087 is a special case. In this 
galaxy the stellar halo starts to affect the shape of the 
surface brightness profile at higher surface brightness lev- 
els. Our results imply that NGC 1087 has an unusually 
bright stellar halo with (^)e ~ 25.5 mag arcsec"^. Prob- 
ably this is the reason why this galaxy , based on shallow 
surface brightness profiles (see IPT06 D , was classified as 
Type III. This altogether raises an important question: 
just how many (late-type) galaxies classified as Type III 
are in reality Type I galaxies with a bright stellar halo? 

The mean effective surface brightness ((/i)e) of our stel- 
lar halos have a typical value of ~ 28 mag arcsec"^ (see 
Table El). NGC 1087 and NGC 7716, the two interacting 
galaxies have brighter (/i)e (see Table [5]). 

An important aspect of our study is the light frac- 
tion contained within the stellar halo compared to the 
total stellar light of the galaxy. We find fractions in 
the range of 1-5%. These are typical values observed 
in nearby galaxies, as well. In two galaxies (NGC 1087 
and NGC 7716), though, the light fraction of the stel- 
lar halo reaches nearly ~ 10%. It is interesting to note 
that in both galaxies we note signs of ongoing interac- 
tion or merging. NGC 7716 is a clear-cut example of 
tidal streams surrounding a galaxy (see Figure [T]), that 
appears as a sharp feature on the surface brightness pro- 
file, and has a surface brightness of ^ 26.5 mag arcsec"^. 

There could be several candidates as progenitors for 
this tidal stream(s), but due to lack of spectroscopic data 
we cannot decide whether the objects along the stream 
do indeed belong to the system or are only background 
objects. There is only o ne known satellite of this galaxy 
reported (Zarit sky et al . 1993), but that is found too far 
(at a projected distance of ~ 130 kpc) to be connected 
to the streams. NGC 1087 is not such an obvious case 
of interaction, but a misaligned outer disk might indi- 



cate that the equilibrium of the disk has been disturbed 
recently by an external source. 

5.2. Color profiles 

The u',g',r',i',and z' filters of SDSS give a full coverage 
of the Spectral Energy Distribution (SED) in the opti- 
cal regime. In principle, using all the five filters gives 
us the opportunity to learn whether the surface bright- 
ness distribution and color gradients exhibit any kind 
of dependence on the observing wavelength allowing us 
to constrain the stellar populations responsible for the 
energy output in the different galactic components. Un- 
fortunately, not all bands behave as well as g' and r'. 
Stacking u'-, i'-, and z'-band data did not result in such 
exceptional images and we cannot derive more than some 
exploratory conclusions. Due to this, we mainly focused 
on g'- and r'-band data. 

We can observe several interesting behavior on the 
published surface brightness and color profiles. Al- 
though, on this sample it is not possible to derive sta- 
tistically robust correlations, we find that the behavior 
of the surface brightness and color is connected to the 
surface brightness types of these galaxies. 

In the Type II-CT galaxies in this sample (like 
NGC 0450 and NGC 0941) we find that the Type II 
behavior appears more pronounced in the u'-band. In 
NGC 0450 the u'-band surface brightness profile seems 
to have a sharp cut-off (see Figure [8]), causing an en- 
hanced "U-shape" on the (u'-g') color profile. Towards 
redder wavelengths the break feature seem to weaken and 
the "U-shape" flattens out (see Figures [51^ 

In BTP08 we obtained robust color profiles of Type I, 
Type II and Type III galaxies. For individual galaxies 
the (g'-r')color profiles reached a depth coinciding with ^ 
25 mag arcsec"^ on the r'-band surface brightness pro- 
files. In the particular case of Type II galaxies, when 
we combined the information of many objects, we were 
able to derive a mean color profile down to ~ 27 mag 
arcsec^^. This mean color profile was not just deeper 
but predicted the general behavior of Type II galaxies. 
Due to that analysis we found that this type of galaxies 
exhibit a "U-shape" on their (g'-r') color. 

Now, using Stripe82 data we got the opportunity to 
trace the color profiles down to fainter levels and confirm 
for instance the classical "U-shape" in case of individual 
galaxies. NGC 0450 and NGC 0941 are typical exam- 
ples of this Type II behavior. Their g'-r' color profiles 
(see Figures [SJ ^ follow nicely the "U-shape" , although, 
the gradient in the inner disk in NGC 0941 is much 
milder but with a prominent reddening beyond the break. 
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The same striking reddening can be observed beyond the 
break in NGC 7716 and UGC 02311. UGC 02081, despite 
of being clasified as Type II based on its surface bright- 
ness profile, exhibits no reddening beyond the break. Ac- 
tuaUy, its flat color profile resambles that of what we 
found to be intrinsic in Type I galaxies. The color pro- 
file of NGC 1068 (Type I) is quite fiat, it has some minor 
reddening around the turn-off radius presumably caused 
by the stellar halo light contribution. The (r'-i') color 
profile of NGC1068 does not reach the same depth, still, 
around the turn-off radius it shows the same slight red- 
dening. 

The color profile of NGC 1087 (see Figure HI]) is puz- 
zling. Based on its surface brightness profile we classified 
this galaxy as Type II*, associating the apparent break 
in the inner disk to the presence of a more concentrated 
star forming region. Meanwhile, the g'-r' color profile ex- 
hibits a fairly deep "U-shape" around that break followed 
with a flattening around the radius where the stellar halo 
contribution starts to domi nate the light. The overall 
shape of this color profile in IBTP08I was associated with 
Type III galaxies. 

Due to the low number of galaxies 10) we did not 
obtain robust color profiles for Type I galaxies in lBTPOH 
although, we have found the same fiattened behavior for 
the mean color profile as along the disk in NGC 1068, 
the only Type I disk we could study in detail in this 
work. NGC 1068 has a rather dusty disk, so it is not 
unexpected that the mean color of the disk of NGC 1068 
is ra ther red (g'-r' ~ 0.6). 

In IBTP08I we found a strong correlation for Type II 
galaxies between the color at the break with the absolute 
magnitude of the galaxy. If we compared the color at the 
break (see in Table |4]) to the the observed (g'-r') = 0.47 
± 0.02 mag mean color of i^TPOl, we should expect that 
the main reason why the majorit y of the o bserved colors 
at the break are bluer than the IBTP08I mean color is 
that these systems bear with less total stellar mass. It is 
indeed true that lower-mass systems such as NGC 0941 
and UGC 02081 exhibit the bluest colors at the break in 
this sample. UGC 02311 also exhibits a blue color at the 
break, despite its high stellar mass. Although, in case of 
this Type II-ORL break one could expect enhanced star 
formation around the resonance zone. The presence of 
star forming rings around the outer Lin dblad Resonance 
has been observed in many s pirals (e.g. iButa fc CrockeJ 
[1991; Buta 2002; G rouchv et al. 201(i). The two pecuhar 
Type II galaxies, NGC 1087 and NG C 7716, a re simflar 
in stellar mass to the mean value of IBTP08I . but only 
NGC 7716 has (g'-r') color close to the BTP08 mean 
value. 

5.3. Stellar surface mass density 

As explained in Section [4. 2. H we obtained M/L ratio 
profiles from (g'-r') colors. These profiles were converted 
into stellar surface mass density (E*) profiles that serve 
as a proxy of the stellar mass distribution of our 7 galax- 
ies. Though, when reaching th e stellar halo the Bell & 
de Jong conversion formula ( Bell fc de ■Tong|[200l seems 
to break down, and the density basically becomes inde- 
pendent of the radial distance. This phenomena is quite 
evident in case of NGC 0450 and NGC 0941 (see e.g. 
Figures [H |9]) . It is important to note that in these two 
galaxies we have observed extremely red colors in their 



stellar halo. If these red colors are due to some extreme 
stel lar populations (with a different IMF, for instance), 
the lBell fc de Jon3 (|2001[ ) empirical formula cannot give 
valid predictions for the M/L, since it is largely based 
on the assumption of a Universal Salpeter-like IMF. 

The profiles of Type II galaxies do not follow the 
shape of the r'-band s urface b rightness profiles. Simi- 
larly to the findings of iBTPOSl the break diminishes on 
the stellar surface mass density profile. To quantify the 
'sharpness' of the break on the profile, we applied 
the same Bulge/Disk decomposition on the E^ profiles 
as for the surface brightness profiles. Not surprisingly, 
the inner-outer scale-lengths ratio takes a value of ~ 1 in 
NGC 0450, NGC 0941 and UGC 02311. In UGC 02081 
the break, however, remains quite prominent, implying a 
different phenomena in this galaxy. (See last column on 
Figure [H) . 

In BTP08 we found that at Type II galaxies, the mean 
stellar mass density at the break radius is ~ 10 Mq 
pc~^. Taking into account that at the break positions 
of NGC 1087 and NGC 7716 the stellar surface mass 
density is much higher than this value (~99 Mq pc~^ 
and ~31 Mq pc~^, respectively), not classifying these 
galaxies as Type II-CT seems justified. On the E^, pro- 
file of UGC 02081 we observe quite an irregular cut-off, 
that seems to be a real mass drop after the break radius. 
Also, the E^, has a low value at the break radius (E^^^br — 
1.7 Mq pc~^), this val ue is much lower than the average 
value found in lBTP08i 

5.4. Stellar populations in the stellar halos 

The observed (g'-r') colors for our stellar halos cover a 
large range from ^ 0.4 to ^ 1.2. This is a rather unex- 
pected variety of colors and stellar populations contained 
within the stellar halos, independently of the total stellar 
mass of the host galaxy. 

In order to learn about the nature of the stellar popu- 
lations behind these observed colors, wc compared them 
to st ellar population model predictions ( Vazdekis et al.l 
l2010f l. As it is already explained in Section l473l Figure 
[4] is only a crude way to visualize what the stellar popu- 
lation content of these stellar halos might be. We chose 
two typical stellar halos to do this analysis, the red stellar 
halo of NGC 1068 ((.g' - r') - 0.55) and the extremely 
red stellar halo of NGC 0941 {{g' - r') ~ 0.95). It is 
important to note, that the observed stellar halo colors 
are redder than that of the disk or even the bulge (with 
the exception of NGC 1068, which has a very dusty disk, 
hence the stellar halo appears to be bluer than the disk). 

It is relatively obvious that with the error bars of 
our measurements neither metallicity, nor age can re- 
ally be constrained. However, by assuming that the 
stellar halo properties of our galaxies are similar to 
the ones in our neighborhood, we can make some fur- 
ther but tentative constraints based on observations of 
nearby galaxies. For instance, it is widely accepted that 
the Milky Way stellar halo is dominat ed by a metal- 
poor population (\Fe/H] < —1.5) (e.g . iRvan fc NorrisI 
Il99lt ICarollo et all \mmn^Je^l6^ and ref- 
erences therein) . In other nearby galaxies, like M31 

llFerguson et al.(l2002l: |Chapmia.n et al]|2006l: llbata et alj 
I 2007 D,^33 CDavidge' '2003': 'McCon nachie et al.l 120061: 
iBarker et al.. .2007. : .Grossi et al.. .20111), or NGC 2403 
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Figure 7. Le/t vniage: NGC7716 in SDSS r'-band (single exposure ) . Ttie limiting surface brightness of this image is ~27 mag aresec" , 
which corresponds to 3(t of the skynoise. The skynoise (similarly to IPTOd ) is the resistant mean of the skyvalues measured on the pixels 
in an aperture surrounding the galaxy. Right image: Stripe82 coaddition of NGC7716 with clearly visible tidal streams likely originating 
from two satellites. The surface brightness of the streams is ^27 mag arcsec"^ in r'-band. 



pavidgell2003HBarker et al.ll2012[ ). the stellar halos con- 
sist of predominantly old and metal-poor {[Fe/H] ^ - 
0.7 - -1.5) stellar populations. 

It is interesting to note that SSP models produced with 
conventional Kroupa IMF (a — 1.3) predict that the stel- 
lar halo color observed in NGC 0941 would imply a stel- 
lar halo that is fairly metal-rich and old. This is at odds 
with the measurements in nearby galaxies. By changing 
the IMF, though, it is possible to look for solutions that 
would attribute the red colors to metal-poor stellar pop- 
ulati ons. It has already be en discussed in the literature 
(e.g. iZackrisson et al.l [20061 and references therein) that 
a bottom-heavy IMF could be responsible for such red 
colors. Following this idea, we applied a bottom-heavy 
IMF with a slope of a = 3.3 to compute our model (g'-r') 
colors. These models indeed yield possible solutions of 
T > 4 Gyr old and [Fe/H] < -1.5. In case of NGC 1068, 
this IMF would imply a young stellar halo (younger than 
in case of Kroupa IMF) with any metallicity possible. 

6. DISCUSSION 

Using deep SDSS Stripe82 data we observed the stellar 
content of 7 late-type galaxies down to Hr' ^ 30 mag 
arcsec"^. This is roughly a stellar surface mass density 
of ~ 0.1 M0 pc~^, allowing us to study the stellar halo 
component of these galaxies. 

At present. Stellar halos are best studied by resolved 
star technique, and there have been only a few at- 
tempts t o study stellar halos by integrated photome- 
try (e.g. iJablonka et al.l |2010[) . However, to explore 
galaxies beyond our neighborhood, integrated photom- 
etry is needed. Resolved star technique provides am- 
ple information on stellar halos: the surface brightness 
of stellar halos in external galaxies, and how their con- 
tribution affects the surface brightness profile. From 
recent observations of nearby galaxies such as M31 or 
NGC2403 we know the structure of the disks and stel- 
lar halos of these galaxies in unprecedented detail (e.e. 



iCourteau eraI1[20Tll: [Barker et al.|[2qT2l ). 

Tha galaxies in our sample are simillar in stellar mass 
and in Hubble-type to many of these well studied sys- 
tems. Although, stochastic variations in the evolutionary 
path of these individual galaxies could cause striking dif- 
ferences as to how the structural properties of the differ- 
ent substructures vary from galaxy to ga laxy on the same 
mass scale. Still, ICourteau et al.l (|2011[ ) argued that one 
should expect to encounter the signatures of stellar ha- 
los in surveys that reach below the de pth of n^i ~27.5 
mag arcsec"^. Using SDSS data lPT06l extracted surface 
brightness profiles down to a similar depth, but the ques- 
tioned remained: What is the real extent of disks ? Does 
the exponential disk continue down to lower surface den- 
sities as predicted by some models? In many cases the 
structure of the outer disk does not change right below 
27 mag arcsec"^. The outer disk indeed follows an expo- 
nential surface brightness profile down to the level where 
there is a smooth continuation into the region where the 
stellar halo light dominates over the disk light (causing 
the upturns discussed beforehand at the level of /i^' ^ 28 
mag arcsec^^). 

Due to this mixture of components, we are unable to 
determine the real extent of the disks. It would seem that 
this disk-stellar halo light conspiracy is due to the observ- 
ing geometry. In face-on geometry we cannot decouple 
the light of the stellar halo from that of the disk. Be- 
cause when we study the surface brightness distribution 
of these two components, we rely on the superimposed 
two dimensional projection of both the stellar halo and 
the disk. The only way to explore further the disk light 
would be in edge-on projection. In edge-on projection by 
integrating the light of the galaxy along the line-of-sight 
of the disk, the contrast between the disk and the stellar 
halo is enhanced, allowing us to explore around ^2 mag 
deeper the disk component. 

Our deep surface brightness profiles has modified our 
view on how surface brightness profiles should be classi- 



16 



fied. By making use of, for instance, high quahty color 
profiles it was possible to rethink the previous classifica- 
tion scheme, especially in relation to Type III profiles. 

6.1. On the nature of Type III profiles 

The nature of T ype III profiles since their discovery 
(|Erwin et al.ll2005l ) proved to be contradictory. Any com- 
prehensive picture regarding the nature of Type III pro- 
files needs to explain: why we observe two sub-class of 
Type III profiles (with disk like or spheroid like isophotes 
in the outskirts), the origin of the peculiar color profile 
of Type III galaxies, and why the statistics of Type III 
profiles decreases towards later Hubble-types. 

The fraction of Type I II profiles becomes larger in 
early- type spirals (|PT06l . lErwin et all 120081 ). Conse- 
quently, as the prominence of the bulge increases with 
the same trend, one could hazard that Type III pro- 
files are nothing else but the bulge component emerg- 
ing above the disk component on the surface brightness 
profiles. This would also explain the roundish isophotes 
that many of these Type III galaxies present in their 
outer regions . Ind eed, being motivated by this idea, 
IMaltbv et al.l (|2012[ ) studied the connection between the 
bulges and Type III profiles in a sample of intermediate 
redshift galaxies. They found that especially in early- 
type disk galaxies Type III profiles are caused by the 
aforementioned bulge emergence, but they could not re- 
fute the idea that in many cases, however, it remains a 
pure disk phenomenon. What is the origin of Type III 
profiles in late-type galaxies where the bulge component 
is rather small to be held responsible for such phenom- 
ena, and where the light of the outer region clearly follows 
disky isphotes? 

In our sample, we have found that the contribution of 
the stellar halo light on the surface brightness profiles 
causes a similar upturn as on 'classical' Type III pro- 
files. In case of NGC 1087 this contribution is so high 
that it starts to affect the surface brightness profile at a 
higher surface brightness level. Due to this, this galaxy 
previously was classified as Type III (PT06; BTP08). It 
is natural to ask whether it is a unique behavior of this 
galaxy or a general rule that would mean that we should 
find many examples once looked carefully. Interestingly, 
M81 (a twin of NGC 1087 with its absolute magnitude 
of Mb — -20.7 mag) is also a s ystem with an unus ually 
bright halo component, though Bar ker et al.l (|2009f ) also 
consider the possibility that the extended bright compo- 
nent in M81 belongs to a disk feature (e.g. thick disk). 
Regardless of this, the light fraction of this component 
compared to the total luminosity in both galaxies is ^ 
10 %. In NGC 7716 the steUar halo is similarly bright, 
moreover in this galaxy the misaligned outer disk, that 
causes a sharp Type III break, can be associated with 
clear interaction. The tidal stream we identified around 
this galaxy is clearly connected to this misaligned outer 
disk. 

NGC 1087 and NGC 7716 have a central star form- 
ing region in their inner disk that causes a high surface 
brightness break on their surface brightness profile that 
we indicated with an "asterix Type 11" profile (distin- 
guishing from the usual Type II classification because 
this break is not the phenomena we usually associate 
with the end of the spiral arms). We can assume that 
the enhanced star formation in the inner disk is trig- 



gered by the recent interaction that has either supplied 
the material necessary to form stars or redistributed the 
gas already settled inside the disk. If this scenario is 
plausible, we can inmediately understand why we ob- 
serve such a peculiar shape of the Type III color profiles. 
As it was already shown in BTPOSi (see their Figure 1) 
Type III galaxies seem to retain the "U-shape" at shorter 
galactocentric radii, with a red "plateau" residing at the 
Type III break followed by a blueing towards the out- 
skirts. It is remarkable that NGC 1087 retains this shape 
on its (g'-r') color profile, at the Type II* break with 
the bluest color and at the Type III break with a red 
"plateau". NGC 7716 also exhibits the bluest color at 
the position of the Type II* break. 

Combining all the information from the surface bright- 
ness and color profiles, we propose the following unified 
picture of Type III (late-type) galaxies: Type III pro- 
files are not an individual class of galaxies but rather. 
Type I or Type II galaxies with a stellar halo. The sur- 
face brightness level at which the stellar halo contributes 
depends on the recent evolutionary history (and stellar 
mass) of the galaxy. Galaxies are transformed by (re- 
cent) interaction. The infalling new material enriches 
the stellar halo, can add extra material to the disk and 
trigger star formation leading to a Type II-I-III mixed 
profile. Meanwhile, a more quiescent history means a 
less developed stellar halo, that would explain the exis- 
te nce of pure Type III profiles ( see the published profiles 
of [PT06l and lErwin et al.l[20?)8h . 

Due to this, in any flux limited survey that cannot 
probe the surface brightness levels of quescient stellar 
halos, the "Type III profiles" would remain unobserved. 
This also woul d expla in the low number of Type II-I-III 
profiles in the IPT06I sample (~ 5%), that necessarily 
comes from the combined effect of a high surface bright- 
ness Type II break and a bright stellar halo/tidal stream 
on the surface brightness profile. Type Il-f-III late-type 
galaxies would be late disks 'caught' during a merger 
event. 

6.2. Further clues on outer-disk formation 

What clues can we deduce from the results of this pa- 
per in relation to the formation of breaks in the surface 
brightness profile of disk galaxies? Can these deep sur- 
face brightness, color and stellar surface mass density 
profiles constrain the formation of these individual sys- 
tems by making use of the predictions of current models? 

The formation of Type II breaks is an issue that has 
been addressed for a long time and with a renewed vigor 
during the past couple of years. As we have already re- 
ferred in Section [1] the new rendition of models incorpo- 
rate the effect of a star formation threshold and the redis- 
tribution of stars by secular processes (migration) as key 
ingredient to match the obser v ational evidence in disk 
galaxies (iRoskar et al.l r2008t ISdnchez-Blazguez et al.l 
120091: iMartmez-Serrano et al.ll2009[ ). These models have 
remarkably similar predictions of an existing minimum 
on the age profile occuring right at the break position 
followed by a likely ageing towards the outskirts. This 
was the natural consequence of the net migration of stars. 

The ageing in the outskirts in case of an absent metal- 
licity gradient can be seen as reddening on the color pro- 
files. The Type II "U-shape" (g'-r') color profiles for this 
reason are in qualitative agreement with the predictions 
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of iRoskar et all (|2008D| iSanchez-Blazguez et all (|2009[ ): 
iMartmez-Serrano et al.l ( 20091 ). A disagree ment betw een 
the observa tional and theor etical res ults of ' BTPOSi and 
IRoskar et al. (2008) was that IBTPOSI proposed the (near) 
absence of the break on the stellar surface mass density 
profile; 

In IBTPOSI we could not show this phenomena on the 
individual galaxy profiles due to the insufficient depth 
of the SDSS data that prevented us to trace the (g'-r') 
color or the S^, profiles well beyond the break rad ius. Al- 
though, based on the robust mean profiles, in IBTPOSI 
we claimed that the break on the surface brightness pro- 
file signifies a radial change in the key ingreadients of the 
stellar population rather than being caused by a mass 
drop beyond the break radius. Based on the results re- 
vealed in this paper, we find strong indication that there 
is no mass drop beyond the break radius. We demon- 
strated on Figure [6] that the inner-out disk scale- lengths 
ratio obtained from the stellar surface mass density drops 
to - 1. 

The inner-outer scale-lengths ratio strongly depends 
on the observing wavelength, as well (see Figure [6]). 
The prominence of the break decreases from u'-band to- 
wards redder filters, in a way that surface brightness 
pro file of the inner disk i s the flattest in the u'-band 
(in IBakos fc Truiillol [20lol we observed the same trend 
on a few selected 'grand-design' spirals.) Consequently, 
the inner disk must be dominated by recently formed, 
young stars that are distributed homogeneously along 
the inner disk, meanwhile old stars can be found ev- 
erywhere in the disk. Remarkably, our observations 
are in qualitative agreement with the predictions of 
ISanchez-Blazguez et al.l (|2009f ) who also showed that the 
break becomes more pronounced in bluer wavebands. It 
is however interesting to note, that they find that migra- 
tion is not the main mechanism producing this particular 
shape of the age profile. The same trend is observed with 
and without migration. When looking for explanation for 
the absence of breaks in the stellar surface mass density 
profile, however, stellar migration is a key process that 
weakens the intensity of the break in the profile. In 
situ star formation in the outer disk cannot produce suf- 
ficient amount of stars to conform the lack of or weak 
breaks. Based on broad-band colors it seems that migra- 
tion and in situ star formation are both relevant. 

On the other hand, we need to stress the importance of 
the discovery of galaxies like UGC 02081 that based on 
their surface brightness profile are classified as Type II 
but do not exhibit the "U-shape" , moreover retain a sig- 
nificant mass drop right after the break radius. By find- 
ing that the E^, at the break is at much lower stellar sur- 
face mass density (~ 1Mqpc~'^) than the observed mean 
value (^ 10Mqpc~^), we think we need to treat them as 
distinct sub-class of Type II profiles. The low E^, value 
at the break and the sharp behavior of the E^, profile of 
UGC 02081 at its outskirts resembles more of the tradi- 
tional picture of truncations (truncations could be imag- 
ined a s the 'edge' of disks) proposed by Ivan der KruitI 
(fl979t) . UGC 02081 could be one of the few moderately 
inclined galaxies where a real truncation is observed. 

Finally, in this paper, we propose that Type III profiles 
are in fact signatures of the stellar halo (or tidal streams) . 
For this reason, the formation of Type III breaks seems 



to be a natural consequence of the current cosmologi- 
cal paradigm. Satellites are expected to interact with 
their host galaxy. Their orbital parameters (angular mo- 
mentum, for instance) will decide how their interaction 
changes the morphology of the host galaxy, whether the 
remnants build the stellar halo or rather settle in the disk 
(jYounger et al.ll2007ll2008t iKazantzidis et al.ll2008l ). We 
iind that in NGC 7716 a disky Type III morphology can 
be associated to a misaligned outer disk that is in connec- 
tion with a tidal stream of an unknown progenitor. The 
jet-like structure more or less perpendicular to the disk 
could be part of another loop. These tidal features are 
well-know a r ound several late-type disks. Just recently, 
iWang et al.l (|2012D investigated the possibility that loops 
in NGC 5907 are formed after a major merger event, chal- 
lenging the dominance of minor merger scenarios to form 
such substructures. 

6.3. Stellar halos beyond the local volume 

We observed the presence of stellar halos in 7 spiral 
galaxies. The signature of the stellar halo is a well- 
defined upturn followed by a flattening on the deep r'- 
band surface brightness profiles, that can be observed in 
the other two deep, g'- and i'-bands. 

Although, the number of systems we can compare to 
is limited, we can find solid examples of stellar halos ob- 
served by 'resolved star' technique in nearby galaxies. If 
we compare the surface brightness profiles of our galaxies 
to star-count profiles of Local Group and other nearby 
galaxies within the next ~15 Mpc, we find that those pro- 
files exhibit the same behavior at similar distances and 
surface brightness levels. In our galaxies the stellar halo 
starts to dominate as far as ~ 10-40 kpc, farther in m ore 
luminous hosts (see Table [5]) . iCourteau et al.l (|2011| ) re- 
ported that in M31 the stellar halo starts to dominate at 
^9 kpc over the light of the disk component. NGC2403 
(Mb=-19.44 mag), a far more isolated system than M31 
and analogous to NGC 0941, also exhibit an extended 
light component that becomes the dominant fea ture of 
the s tellar count profile at ~18 kpc radius (Barke r et al.l 
[20T1 ). 

The light fraction of the stellar halos in our galaxies 
fall in the range of 1 - 5 %. These are typical values 
found in nearby galaxies, the stellar halo light fractions 
observed in the Milky Way (jCarollo et a l. 2010), M31 
(Cou rteau etall l201llf ) or in NGC 2403 (Barker eLaU 
2012) are of the same order we find in our galaxies. The 
light fraction found in an asymmetric substructu re sur- 
rounding M33 is ~ 1% (iMcConnachie et al.llMol) 

According to the findings of iPurcell et al.l ([20071 ) 
that investigated the formation of stellar halos through 
merger-history, the light fraction in stellar halos depends 
on the total stellar mass of the system and the accre- 
tion history of the galaxy, hence on the mass of the dark 
matter halo. Although, the observed stellar mass (or ab- 
solute magnitude) range of our galaxies is of the order 
of 4 magnitudes, we do not find evidence of the stellar 
halo light increasing with the host galaxy mass. This 
might occur if due to the relatively low number of galax- 
ies observed in this sample, the variations in the accre- 
tion history of these galaxies smears out the correlation 
expected to be seen with the total galaxy mass. We do 
find that in the two systems where we discovered the 
signs of recent interaction the stellar halo light fraction 
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is systematically higher 10%). 

Other cosmological simulations also emphasize the pos- 
sibility of a dual origin of stellar halos. In this scenario, 
stellar halos should contain stars that were formed in 
situ and stars that were accrete d from lower ma ss satel- 
lites (e.g. .MotmLSialj .^^eEiD Iloiil). Con- 
sequently, more massive galaxies are expected to have 
more metal rich stellar halos. On the o bservational side 
there are contradictory findings so far. iMouhcine et al.l 
suggested that the metallicity of stellar halos 
in more massive sys tems is indeed higher, however, 
iFerguson et al.l (|2007[ ) showed that the observed photo- 
metric metallicities can be largely biased that can be 
overcome by kinematically selecting the halo star popu- 
lation. As a results of such a careful treatment of halo 
stars, they found no evidence of correlation between stel- 
lar halo metallicity and host galaxy luminosity. 

The stellar halo colors we find in our galaxies show no 
correlation with the host galaxy mass (absolute magni- 
tude) either. We observed a range of different halo col- 
ors from red to extremely red, but we have found both 
'blue' and red colors in lower-mass systems, namely in 
NGC 0941 (See Figure [9]) and UGC 02081 (See Figure 
[9]) . If not the host galaxy mass then what drives the color 
distribution of these stellar halos? Do galaxies evolving 
in isolation have redder colors? Is there any reason to be- 
lieve that a variant IMF is responsible for the observed 
colors? Or is there a large population of stars formed in 
situ that could account for the presence of an extremely 
old but metal-poor population? These are questions that 
so far remain without answer. 

We used the latest MILES SSP models to give crude 
constraints on the stellar populations embedded within 
these stellar halos. It is evident that resolved star tech- 
nique yields more accurate stellar population properties, 
like metallicities. Broad-brand colors are affected by the 
age-metallicity degeneracy, so it is not possible to ascer- 
tain the stellar populations dominating the stellar halos 
of our galaxies. The typical stellar halo colors seem be 
in agreement with the ages and metallicities of nearby 
galaxies (r ~ 5 Gyr and [Fe/H] ^ -1.5). But extreme red 
halo colors cannot be reconciled with such populations or 
with the current cosmological paradigm. We find com- 
pelling evidence that these extreme red colors may indi- 
cate a varying IMF. A bottom-heavy IMF could be the 
key agent in producing such red populations. Such solu- 
tion to the red halo phe nomenon has also been sugge sted 
by several authors (e.g. IZackrisson et al.ll2006l 120121 ). 

7. CONCLUSIONS 

1. The light contribution from the stellar halo in the 
surface brightness profiles of spiral galaxies causes 
Type III profiles: the surface brightness level at 
which the surface brightness profile starts to devi- 
ate from the exponential profile of the disk depends 
on the brightness of the stellar halo. The brighter 
the stellar halo, the brighter the Type III "break" . 
Disky Type III profiles are caused by tidal streams, 
that also belong to the stellar halo. These effects 
could be responsible for previous classification of 
surface brightness profiles of late-type galaxies as 
Type III. 

2. High surface brightness disk breaks {fir' ~22 mag 



arcsec"^) could be an indication of recent interac- 
tion. The new infalling material (gas) settles into 
the center of the disk and triggers star formation 
in the inner disk. 

3. The stellar halos of quiescent systems contribute 
less to the total light of the galaxy. 

4. Stellar halos show a large range of {g' — r') colors. 
Some of them, are extremely red {g' — r' >1), al- 
lowing for the possibility of bottom-heavy IMFs to 
explain the red halo colors. 

8. SUMMARY 

In this paper we have explored the structure and stel- 
lar population content of 7 late-type spirals using re- 
stacked deep data from SDSS Stripe82. This data typi- 
cally reaches ~ 2 magnitudes deeper in the 5 bands than 
regular SDSS data. In the r'-band the quality of the data 
made it possible to reach a surface brightness level of ^ 
30 mag arcsec"^. 

Using the deep r'-band profiles we classified our galax- 
ies into the three classical types (Type I, Type II, and 
Type III). All surface brightness profiles change their be- 
havior around ^ 28 mag arcsec"^. A flattening is ob- 
served on all profiles that we believe to be the signature 
of stellar halos. 

These faint levels composed of the outer stellar disk 
and stellar halo beforehand were only reached by tech- 
niques that involve resolving stars in these regions or 
stacking many individual objects. The resolved star tech- 
nique, although is very powerful, is limited to galaxies 
that are in the reach of the resolving power of current 
telescopes. Studies of nearby galaxies like M31, M33, 
NGC2403 showed that the faint outer regions consist of 
an ample stellar content with many substructures present 
in them. We have found different stellar halo color im- 
plying very distinct stellar halo populations among our 
galaxies. We find compelling evidence that there exist 
extreme red colors halo colors that can be explained by 
stellar populations that are formed from a bottom-heavy 
IMF. Meanwhile, less red colors we find in other stellar 
halos could imply metal-poor ([Fe/H] ~ -1.5) and old 
(r > 5 Gyr) stellar populations consistent with the stel- 
lar populations found in more nearby galaxies. The light 
fraction in the stellar halos of our galaxies is 1-5%, but 
in galaxies where recent interaction occured, it can reach 
10%. In NGC 7716 we have discovered the presence of a 
giant stellar stream, and a ~ -14 r'-mag satellite in the 
vicinity of NGC 1068. 
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APPENDIX 
A. GALAXY ATLAS 

In this section we publish the surface brightness, color, stellar surface mass density profiles of the 7 spiral galaxies 
of our sample. These figures contain a RGB color composite stamp created from the deep g'-,r'-, and i'-band stacks. 
These stacks were created using STIFF. In order to help the identification of low surface brightness features by naked 
eye, we inverted the black background later on. In some cases we also used "unsharp masking" to enhance these 
features. (All images have the same orientation with North displayed upwards and with East to left.) 

The top right panel contains the r'-band profile of our galaxies. The plot contains the observed profiles, and the 
lcr(blue shade) and 3cr(purple shade) error regions. We drew a blue dashed line to indicate the 3cr noise-level of the 
sky, the radius where the profile reaches this level is indicated with the blue ellipse on the color image. 

The r'-band profile served as the basis for decomposing the different structures building up the galaxies, like bulge 
{green dots), disk (black dotted- dashed lines), and stellar halo (yellow dots) components. The decomposition is ex- 
plained in detail in Section 14.1.31 The final composite profile is drawn over the observed profile by red open circles. 
The region where the stellar halo contributes to the total light of the galaxy with more than 50% is indicated by the 
green shade. If the surface brightness profile of the galaxy exhibits a break, it is shown by a vertical red dashed line 
on the profile, and with a red ellipse on the color image. 

The middle left hand panel shows the (g'-r') color profile using scarcely placed bins. For the stellar halo we derived 
one single color. 

The middle right hand panel shows the S^, profile (for further information see Section l4.2.3p . 

The bottom right and left panels contain the surface brightness profiles and color profiles in all available SDSS filters 
obtained from the Stripe82 imaging. The surface brightness and color profiles are color coded as indicated on the 
figures. 
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Fi gure 8 . NGC 0450: It is a moderately inclined spiral galaxy with a background galaxy (UGC00807) superimposed. It has been classified 
bv VPTOg as a Type II-CT. 
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Figure 9. NGC 0941: It is a moderately inclined SABc type galaxy. It is within a small group of galaxies, though seems remarkably 
isolated and the disk shows no signs of interaction. It has a small bar that generates the pseudobulge feature observed at R < 30" 
(Kormendy, private communication). The break lies at ~ 70 arcseconds (~ 7.5 kpc), and is classified as Type II-CT. 
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Figu re 10. NGC 1068 or M77: This galaxy is a frequently studied, bright galaxy, wit h a Sey f ert 2 c lass AGN IjKhachikian & Weedmar] 
119741) in its center. NGC 1068 appears to have a complex inner structure. According to lErwinI }2004!) this is a double-barred galaxy with 
a small inner bar (R < 20") e mbedded within an oval distortion (extending out to R < 200"; Schinncrcr ct al. 2000) also remarked by 
IKormendv fc KennicuttI II2004I ) . We consider the structure outside R > 200" as the genuine disk of the galaxy. This galaxy is basically 
analogous to NGC 4736 (M94) (see also e.g. ITruiillo et al . 2009) showing a strikingly similar overall morphology. Considering all these 
morphological features, this galaxy, based on its r'-band surface brightness profile, has been classified as Type I. The extended outer disk 
shows signs of interaction and misalignment, probably a result of recent bombardment. We report here the discovery of a low surface 
brightness (/i^/ ~27 mag arcsec"'^) satellite found at the coordinates a = 02 : 43 : 00 5 = —00 : 09 : 01. This object is indicated with a red 
arrow. 
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Figure 11. NGC 1087: It is a bright, SABc type galaxy with a misaligned outer disk that shows signs of an ongoing interaction on the 
North side of the galaxy. NGC 1087 is classified as Type II*+III: the apparent Type II break is due to the central start forming region 
in the inner disk, meanwhile the bright Type III like feature starting at R ~ 120" appears because of the presence of an unusually bright 
stellar halo (due to this, this galaxy previously was classified as Type III in lPTOg ). 
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Figure 12. NGC 7716: One of the most interesting galaxies of our sample. This galaxy shows no pecularities on its DR7 image, but its 
Stripe82 counterpart reveals the presence of a misaligned outer disk and stellar streams. We cannot identify the progenitor satellites of 
these streams from the images (curiously enough the object that looks like a bright satellite to the North of NGC 7716 is just a background 
galaxy). Unfortunately, the nearby brights star obscures a part of the stellar stream. The other, jet-like feature is possibly another stream 
projected along our light-of-sight. NGC 7716 is classified as Type following the previous classification scheme. The Type III 

feature, however, is clearly linked to the presence of the bright stellar stream and is part of the stellar halo. 
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Figure 13. UGC 02081: Faint, moderately inclined SABc spiral. The bar is not obvious from the image, nor from the profile. It has a 
faint break at /i^/ ~ 25 mag arcsec"^ and is classified as Type II-CT. 
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Figure 14. UGC 02311: Bright galaxy with a strong bar and with an apparent star forming resonance ring. The break position coincides 
with the position of the ring, therefore this galaxy is classified as Type II-ORL. 
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B. SURFACE BRIGHTNESS PROFILES 

In this section we publish the surface brightness profiles of our galaxies obtained from the Stripe82 imaging in the 
■5 SDSS bands. For some of our galaxies, the u'-band and z'-band detector problems prevented us to create reliable 
mosaics. In these cases (NGC 1068, NGC 1087, UGC 02081, and UGC 02311 in the u'-band and NGC 1068 in the 
z'-band ) we did not use the u'- and z'-band stacks for further analysis. 



Table 6 

SurfELce brightness profiles of NGC0450 



Radius 
[arcsec] 


[mag arcsec"^] 


[mag arcsec"^] 


r — 2i 

[mag arcsec J 


Mi' 

[mag arcsec"^] 


[mag arcsec J 


2.0 


21 


44 ± 0.02 


20 


26 ± 0.02 


1 Q en + n 01 


19 


55 ± 0.01 


in QQ + n ni 

X£7.(J^ _1_ \J.\JJ. 


5.9 


22 


08 ± 0.02 


20 


90 ± 0.01 


20.40 ±0.01 


20 


15 ± 0.01 


19.92 ±0.01 


10.0 


22 


67 ±0.01 


21 


43 ±0.01 


20.92 ± 0.01 


20 


66 ± 0.01 


20.45 ± 0.01 


14.1 


22 


90 ± 0.01 


21 


77 ±0.01 


21.26 ± 0.01 


21 


03 ±0.01 


20.83 ±0.01 


18.4 


23 


12 ±0.01 


22 


03 ± 0.00 


21.56 ±0.00 


21 


32 ± 0.00 


21.13 ±0.00 


22.8 


23 


23 ± 0.01 


22 


20 ± 0.00 


21.78 ±0.00 


21 


55 ± 0.00 


21.38 ± 0.00 


27.3 


23 


40 ± 0.01 


22 


39 ± 0.01 


21.98 ± 0.01 


21 


77 ± 0.00 


21.61 ±0.01 


31.9 


23 


62 ± 0.01 


22 


60 ± 0.01 


22.18 ± 0.01 


21 


98 ± 0.00 


21.82 ± 0.01 


36.7 


23 


66 ±0.01 


22 


66 ± 0.01 


22.27 ±0.01 


22 


09 ±0.01 


21.93 ±0.01 


41.7 


23 


76 ± 0.01 


22 


73 ± 0.01 


22.36 ±0.01 


22 


19 ±0.01 


22.04 ±0.01 


46.8 


23 


84 ± 0.01 


22 


81 ±0.00 


22.47 ± 0.00 


22 


31 ± 0.00 


22.16 ±0.01 


52.0 


23 


99 ± 0.01 


22 


96 ± 0.00 


22.63 ± 0.00 


22 


47 ± 0.00 


22.33 ± 0.00 


57.4 


24 


22 ± 0.01 


23 


17 ± 0.00 


22.84 ± 0.00 


22 


67 ±0.00 


22.54 ± 0.00 


63.0 


24 


34 ±0.01 


23 


36 ± 0.00 


23.03 ± 0.00 


22 


87 ± 0.00 


22.74 ±0.00 


68.7 


24 


41 ±0.01 


23 


51 ±0.00 


23.22 ±0.00 


23 


07 ±0.00 


22.94 ±0.00 


74.6 


24 


47 ±0.01 


23 


61 ±0.00 


23.34 ± 0.00 


23 


22 ± 0.00 


23.09 ± 0.01 


80.7 


24 


76 ± 0.01 


23 


85 ± 0.00 


23.56 ± 0.00 


23 


44 ± 0.00 


23.32 ±0.01 


87.0 


25 


16 ±0.01 


24 


19 ± 0.01 


23.88 ± 0.00 


23 


76 ± 0.01 


23.66 ± 0.01 


93.4 


25 


59 ± 0.02 


24 


56 ±0.01 


24.20 ±0.01 


24 


06 ± 0.01 


23.93 ± 0.01 


100.1 


25 


95 ± 0.02 


24 


94 ± 0.01 


24.58 ± 0.01 


24 


44 ± 0.01 


24.32 ±0.01 


106.9 


26 


56 ± 0.03 


25 


36 ± 0.01 


24.98 ± 0.01 


24 


86 ± 0.01 


24.72 ± 0.01 


113.9 


27 


46 ± 0.05 


25 


88 ± 0.01 


25.42 ± 0.01 


25 


33 ± 0.01 


25.17 ±0.02 


121.2 


28 


70 ±0.14 


26 


42 ± 0.01 


25.90 ±0.01 


25 


86 ±0.01 


25.64 ±0.03 


128.7 






26 


94 ±0.01 


26.33 ± 0.01 


26 


28 ± 0.02 


26.15 ± 0.04 


136.4 






27 


38 ±0.02 


26.72 ± 0.01 


26 


82 ± 0.02 


26.57 ± 0.05 


144.3 






27 


78 ± 0.03 


27.02 ± 0.02 


27 


27 ± 0.03 


26.89 ± 0.07 


152.5 






28 


16 ± 0.04 


27.22 ± 0.02 


27 


47 ± 0.04 


27.43 ±0.11 


160.9 






28 


51 ± 0.04 


27.53 ±0.02 


27 


99 ± 0.06 




169.6 






28 


86 ±0.06 


27.79 ± 0.03 


28 


30 ± 0.07 




178.5 






29 


20 ± 0.08 


27.92 ± 0.03 








187.7 










28.20 ±0.04 








197.2 










28.41 ±0.04 








206.9 










28.57 ±0.05 








217.0 










28.76 ± 0.06 








227.3 










28.81 ± 0.06 








238.0 










29.10 ±0.08 
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Table 7 

SurfELce brightness profiles of NGC0941 



Radius Hy^i Hgi /i^/ /i^/ 

[arcsec] [mag arcsec"^] [mag axcsec"^] [mag axcsec"^] [mag arcsec"^] [mag arcsec"^] 



2.0 


21 


84 


± 0.02 


20 


83 ± 


02 


20 


33 


± 


01 


20 


08 ± 


01 


19 


95 


± 


01 


5.9 


22 


20 


± 0.01 


21 


12 ± 


01 


20 


62 


±0 


01 


20 


37 ± 


00 


20 


26 


±0 


00 


10.0 


22 


33 


± 0.01 


21 


30 ± 


01 


20 


84 


± 


01 


20 


61 ± 


00 


20 


51 


± 


01 


14.1 


22 


65 


± 0.01 


21 


59 ± 


01 


21 


15 


±0 


01 


20 


93 ± 


01 


20 


85 


±0 


01 


18.4 


23 


03 


±0.01 


21 


95 ± 


01 


21 


51 


±0 


01 


21 


27 ±0 


00 


21 


19 


±0 


00 


22.8 


23 


25 


±0.01 


22 


24 ±0 


00 


21 


81 


±0 


00 


21 


59 ±0 


00 


21 


51 


±0 


00 


27.3 


23 


58 


± 0.01 


22 


51 ± 


00 


22 


10 


± 


00 


21 


88 ±0 


00 


21 


81 


±0 


00 


31.9 


23 


82 


± 0.01 


22 


77 ±0 


00 


22 


38 


± 


00 


22 


17 ± 


00 


22 


11 


± 


00 


36.7 


24 


21 


± 0.01 


23 


16 ± 


00 


22 


75 


±0 


00 


22 


53 ± 


00 


22 


48 


±0 


00 


41.7 


24 


45 


±0.01 


23 


42 ± 


00 


23 


01 


±0 


00 


22 


79 ±0 


00 


22 


73 


±0 


00 


46.8 


24 


70 


± 0.01 


23 


63 ± 


00 


23 


24 


±0 


00 


23 


01 ± 


00 


22 


96 


±0 


00 


52.0 


24 


97 


± 0.01 


23 


92 ± 


00 


23 


53 


± 


00 


23 


31 ± 


00 


23 


26 


± 


01 


57.4 


25 


27 


± 0.01 


24 


22 ± 


00 


23 


84 


±0 


00 


23 


61 ± 


00 


23 


57 


±0 


01 


63.0 


25 


51 


± 0.01 


24 


51 ± 


00 


24 


12 


±0 


00 


23 


89 ±0 


00 


23 


85 


±0 


01 


68.7 


25 


79 


± 0.01 


24 


81 ± 


00 


24 


42 


± 


00 


24 


18 ± 


00 


24 


13 


± 


01 


74.6 


26 


14 


± 0.02 


25 


24 ±0 


01 


24 


81 


±0 


01 


24 


57 ±0 


01 


24 


49 


±0 


01 


80.7 


26 


70 


±0.02 


25 


77 ±0 


01 


25 


31 


±0 


01 


25 


08 ±0 


01 


25 


01 


±0 


02 


87.0 


27 


38 


±0.04 


26 


38 ±0 


01 


25 


86 


±0 


01 


25 


63 ±0 


01 


25 


54 


±0 


03 


93.4 


27 


99 


±0.07 


26 


94 ±0 


01 


26 


39 


±0 


01 


26 


17 ±0 


02 


26 


03 


±0 


04 


100.1 


28 


99 


±0.16 


27 


58 ±0 


02 


26 


95 


±0 


02 


26 


81 ±0 


03 


26 


56 


±0 


06 


106.9 








28 


07 ±0 


03 


27 


43 


±0 


02 


27 


28 ± 


03 


27 


00 


±0 


09 


113.9 








28 


70 ±0 


05 


27 


80 


±0 


03 


27 


91 ±0 


06 


27 


31 


±0 


11 


121.2 








29 


35 ±0 


09 


28 


29 


±0 


05 


28 


44 ±0 


09 










128.7 














28 


75 


±0 


07 


28 


93 ±0 


13 










136.4 














29 


05 


±0 


09 
















144.3 














29 


33 


±0 


10 

















31 



Table 8 

SurfELce brightness profiles of NGC1068 



Radius 




Hg' 












Mi' 




[arcsecl 


[mag axcsec ] 


[mag axcsec ] 


[mag arcsec 


2.0 


16, 


.35 ± 0, 


,03 


15, 


,78 ± 0, 


,02 


15. 


,58 it 0, 


.02 


5.9 


17. 


,49 ± 0, 


,01 


16, 


,71 ± 0, 


,01 


16. 


,40 it 0, 


.01 


10.0 


18, 


,06 ± 0, 


,01 


17, 


,39 ± 0, 


,01 


17. 


,08 it 0, 


.01 


14.1 


18, 


,61 ± 0, 


,01 


17, 


,92 ± 0, 


,01 


17. 


,62 it 0, 


.01 


18.4 


19, 


.36 ± 0, 


.01 


18, 


.56 ± 0, 


,01 


18. 


,23 it 0, 


.01 


22.8 


19, 


,84 ± 0, 


,01 


19, 


,06 ± 0, 


,01 


18, 


,70 it 0, 


.00 


27.3 


20, 


,00 ± 0, 


,00 


19, 


,30 ± 0, 


,00 


18, 


,96 it 


.00 


31.9 


20, 


,08 ± 0, 


,00 


19, 


,44 ± 0, 


,00 


19, 


,13 it 


.00 


36.7 


20, 


,27 ± 0, 


,00 


19, 


,65 ± 0, 


,00 


19, 


,34 ± 


.00 


41.7 


20, 


,53 ± 0, 


,00 


19, 


,90 ± 0, 


,00 


19, 


,58 it 


.00 


46.8 


20, 


,75 ± 0, 


,00 


20, 


,11 ± 0, 


,00 


19, 


,79 it 


.00 


52.0 


21, 


.00 ± 0, 


.00 


20, 


.33 ± 0, 


.00 


20, 


,00 it 0, 


.00 


57.4 


21, 


.26 ± 0, 


.00 


20, 


.57 ± 0, 


.00 


20, 


,22 it 


.00 


63.0 


21, 


.51 ± 0, 


.00 


20, 


.81 ± 0, 


.00 


20, 


,45 it 0, 


.00 


68.7 


21, 


.72 ± 0, 


.00 


21, 


.05 ± 0, 


.00 


20, 


,70 it 0, 


.00 


74.6 


21, 


.94 ± 0, 


.00 


21, 


.30 ± 0, 


.00 


20, 


,97 it 


.00 


80.7 


22, 


,28 ± 0, 


,00 


21, 


,63 ± 0, 


,00 


21, 


,29 it 0, 


.00 


87.0 


22 


,60 ± 


,00 


21 


,95 ± 


,00 


21, 


,60 it 0, 


.00 


93.4 


22 


,86 ± 


,00 


22 


,21 it 


,00 


21, 


,86 it 


.00 


100.1 


23 


,06 ± 


,00 


22 


,41 ± 


,00 


22, 


,06 it 


.00 


106.9 


23 


,18 ± 


,00 


22 


,54 ± 


,00 


22, 


,20 it 


.00 


113.9 


23, 


,27 ± 0, 


,00 


22 


,62 ± 0, 


,00 


22, 


,28 it 


.00 


121.2 


23, 


,32 ± 0, 


.00 


22, 


.68 ± 0, 


.00 


22, 


,35 it 0, 


.00 


128.7 


23, 


.39 ± 0, 


.00 


22, 


.76 ± 0, 


.00 


22, 


,42 it 


.00 


136.4 


23, 


.48 ± 0, 


.00 


22, 


.85 ± 0, 


.00 


22, 


,51 it 0, 


.00 


144.3 


23, 


.57 ± 0, 


.00 


22, 


.94 ± 0, 


.00 


22, 


,60 it 


.00 


152.5 


23, 


.65 ± 0, 


.00 


23, 


.01 ± 0, 


.00 


22, 


,68 it 


.00 


160.9 


23, 


,69 ± 0, 


,00 


23, 


,06 ± 0, 


,00 


22, 


,73 it 0, 


.00 


169.6 


23, 


,72 ± 0, 


,00 


23, 


,11 ± 0, 


,00 


22, 


78 ± 0, 


,00 


178.5 


23, 


,80 ± 0, 


,00 


23, 


,19 ± 0, 


,00 


22, 


,86 ± 0, 


,00 


187.7 


23, 


,94 ± 0, 


,00 


23, 


,33 ± 0, 


,00 


22, 


,99 ± 0, 


,00 


197.2 


24, 


,07 ± 0, 


,00 


23, 


,46 ± 0, 


,00 


23, 


12 ± 0, 


,00 


206.9 


24, 


,24 ± 0, 


,00 


23, 


,63 ± 0, 


,00 


23, 


29 ± 0, 


,00 


217.0 


24, 


,43 ± 0, 


,00 


23, 


,81 it 0, 


,00 


23, 


,46 it 0, 


.00 


227.3 


24, 


,60 ± 0, 


,00 


23, 


,98 it 0, 


,00 


23. 


,62 it 0, 


.00 


238.0 


24, 


,74 ± 0, 


,00 


24, 


,12 it 0, 


,00 


23. 


,75 it 0, 


.00 


249.0 


24, 


,88 ± 0, 


,00 


24, 


,26 it 0, 


,00 


23. 


,89 it 0, 


.00 


260.3 


25, 


,01 ± 0, 


,00 


24, 


,39 it 0, 


,00 


24, 


,01 it 0, 


.00 


271.9 


25, 


,20 ± 0, 


,00 


24, 


,57 ± 0, 


,00 


24, 


,19 it 0, 


.00 


283.9 


25, 


,41 ± 0, 


,00 


24, 


,77 ± 0, 


,00 


24, 


,40 ± 0, 


.00 


296.3 


25, 


,62 ± 0, 


,00 


24, 


,99 ± 0, 


,00 


24, 


60 ± 0, 


,00 


309.0 


25, 


,87 ± 0, 


,00 


25, 


,23 ± 0, 


,00 


24, 


84 ± 


,00 


322.1 


26, 


,09 ± 0, 


,00 


25, 


,46 ± 0, 


,00 


25, 


07 ± 


,00 


335.6 


26, 


,32 ± 0, 


,00 


25, 


,68 ± 0, 


,00 


25, 


29 ± 


,00 


349.5 


26, 


,55 ± 0, 


,00 


25, 


,91 it 0, 


,00 


25, 


,51 it 0, 


.00 


363.9 


26, 


,74 ± 0, 


,00 


26, 


,11 it 0, 


,00 


25. 


,69 it 0, 


.00 


378.6 


26, 


,99 ± 0, 


,00 


26, 


,36 it 0, 


,00 


25. 


,93 it 


.00 


393.8 


27, 


,20 ± 0, 


,00 


26, 


,58 it 0, 


,00 


26, 


,12 it 


.00 


409.5 


27, 


,49 ± 0, 


,00 


26, 


,84 it 0, 


,00 


26, 


,36 it 0, 


.00 


425.6 


27, 


,69 ± 0, 


,01 


27, 


,04 ± 0, 


,00 


26, 


,55 ± 0, 


.00 


442.2 


27, 


,97 ± 0, 


,01 


27, 


,29 ± 0, 


,01 


26, 


75 ± 0, 


.01 


459.3 


28, 


,17 ± 0, 


,01 


27, 


,46 ± 0, 


,01 


26, 


,93 ± 0, 


.01 


476.9 


28, 


,31 ± 0, 


,01 


27, 


,59 ± 0, 


,01 


27, 


,04 ± 0, 


,01 


495.1 


28, 


,44 ± 0, 


,01 


27, 


,71 ±0, 


,01 


27, 


17±0, 


,01 


513.8 


28, 


,53 ±0, 


,01 


27, 


,83 ±0, 


,01 


27, 


,28 ± 0, 


,01 


533.0 


28, 


,50 ±0, 


,01 


27, 


,88 it 0, 


,01 


27. 


,30 it 0, 


.01 


552.9 


28, 


,62 ± 0, 


,01 


28, 


,04 it 0, 


,01 


27. 


,44 it 0, 


.01 


573.3 


28, 


,74 ± 0, 


,01 


28, 


,15it0. 


,01 


27. 


,52 it 0, 


.01 


594.3 


28, 


,89 ± 0, 


,01 


28, 


,39 it 0, 


,01 


27, 


,73 it 0, 


.01 


616.0 


29, 


,05 ±0, 


,01 


28, 


,.58 ± 0, 


,01 


27, 


,94 ±0 


,01 


638.3 


29, 


,18 ±0, 


,01 


28, 


,81 ±0, 


,02 


28, 


12 ± 0, 


,01 


661.3 


29, 


,25 ±0, 


,01 


28, 


,92 ±0, 


,02 


28, 


15 ± 0, 


.01 


685.0 


29, 


,59 ±0, 


,02 


29, 


,29 ±0, 


,02 


28, 


,27 it 0, 


.02 


709.4 


29, 


,82 ±0, 


,02 


29, 


,47 ± 0, 


,03 


28. 


,50 it 0, 


.02 


734.5 


30, 


,64 ± 0, 


,05 


30, 


,71 ±0, 


,08 
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Table 9 

Surface brightness profiles of NGC1087 



Radius 






Mi' Mz' 


[arcsec] 


[mag arcsec 


[mag arcsec"^] 


[mag arcsec"^] [mag arcsec"^] 



2.0 


19 


25 ± 


04 


18 


70 ± 


04 


18 


53 ± 


03 


18 


25 + 03 


5.9 


20 


35 ±0 


01 


19 


71 ± 


01 


19 


37 ±0 


01 


19 


10 it 0.01 


10.0 


20 


60 ± 


01 


20 


02 ± 


01 


19 


73 ± 


01 


19 


50 ± 0.01 


14.1 


20 


70 ± 


01 


20 


18 ± 


01 


19 


94 ± 


00 


19 


74 ± 0.01 


18.4 


20 


74 ±0 


00 


20 


26 ±0 


00 


20 


05 ±0 


00 


19 


87 ± 0.00 


22.8 


20 


71 ±0 


00 


20 


28 ±0 


00 


20 


08 ±0 


00 


19 


92 ± 0.00 


27.3 


20 


91 ± 


01 


20 


51 ± 


00 


20 


29 ±0 


00 


20 


15 ± 0.00 


31.9 


21 


16 ± 


00 


20 


74 ± 


00 


20 


52 ± 


00 


20 


38 ± 0.00 


36.7 


21 


24 ± 


00 


20 


84 ±0 


00 


20 


62 ±0 


00 


20 


48 ± 0.00 


41.7 


21 


37 ±0 


00 


20 


95 ± 


00 


20 


73 ±0 


00 


20 


59 ± 0.00 


46.8 


21 


52 ± 


00 


21 


09 ± 


00 


20 


87 ± 


00 


20 


74 ± 0.00 


52.0 


21 


71 ± 


00 


21 


29 ± 


00 


21 


08 ± 


00 


20 


94 ± 0.00 


57.4 


21 


91 ± 


00 


21 


50 ± 


00 


21 


28 ± 


00 


21 


15 ± 0.00 


63.0 


22 


20 ±0 


00 


21 


78 ±0 


00 


21 


54 ±0 


00 


21 


42 ± 0.00 


68.7 


22 


50 ± 


00 


22 


08 ± 


00 


21 


82 ±0 


00 


21 


72 ± 0.00 


74.6 


22 


78 ± 


00 


22 


34 ± 


00 


22 


08 ± 


00 


21 


98 ± 0.00 


80.7 


23 


01 ± 


00 


22 


57 ±0 


00 


22 


30 ±0 


00 


22 


20 ± 0.00 


87.0 


23 


33 ±0 


00 


22 


87 ±0 


00 


22 


58 ±0 


00 


22 


48 ± 0.00 


93.4 


23 


64 ±0 


00 


23 


16 ±0 


00 


22 


86 ±0 


00 


22 


76 ± 0.00 


100.1 


23 


93 ±0 


00 


23 


42 ±0 


00 


23 


11 ±0 


00 


23 


01 ± 0.00 


106.9 


24 


20 ± 


00 


23 


68 ± 


00 


23 


35 ± 


00 


23 


26 ± 0.00 


113.9 


24 


48 ± 


00 


23 


97 ±0 


00 


23 


63 ±0 


00 


23 


53 ± 0.00 


121.2 


24 


73 ±0 


00 


24 


20 ±0 


00 


23 


87 ±0 


00 


23 


75 ± 0.01 


128.7 


25 


01 ± 


00 


24 


49 ± 


00 


24 


14 ± 


00 


24 


07 ± 0.01 


136.4 


25 


25 ± 


00 


24 


74 ± 


00 


24 


39 ± 


00 


24 


33 ± 0.01 


144.3 


25 


48 ± 


00 


24 


97 ± 


00 


24 


63 ± 


00 


24 


57 ± 0.01 


152.5 


25 


70 ± 


00 


25 


20 ± 


00 


24 


86 ±0 


00 


24 


80 ± 0.01 


160.9 


25 


88 ± 


00 


25 


39 ± 


00 


25 


05 ±0 


00 


25 


00 ± 0.01 


169.6 


26 


09 ± 


00 


25 


60 ± 


00 


25 


28 ± 


01 


25 


23 ± 0.02 


178.5 


26 


35 ± 


01 


25 


85 ± 


01 


25 


55 ± 


01 


25 


47 ± 0.02 


187.7 


26 


58 ±0 


01 


26 


14 ±0 


01 


25 


83 ±0 


01 


25 


82 ± 0.03 


197.2 


26 


90 ±0 


01 


26 


45 ±0 


01 


26 


24 ±0 


01 


26 


26 ± 0.05 


206.9 


27 


20 ±0 


01 


26 


79 ±0 


02 


26 


55 ± 


02 


26 


65 ± 0.07 


217.0 


27 


38 ± 


02 


26 


89 ± 


02 


26 


74 ± 


02 


26 




227.3 


27 


58 ±0 


02 


27 


11±0 


02 


27 


05 ±0 


03 


27 


26 ±0.11 


238.0 


27 


88 ±0 


02 


27 


35 ±0 


02 


27 


35 ±0 


04 


27 


98 ± 0.20 


249.0 


28 


17 ±0 


03 


27 


66 ±0 


03 


27 


57 ±0 


04 


28 


08 ± 0.20 


260.3 


28 


40 ±0 


03 


27 


87 ±0 


03 


28 


05 ±0 


06 






271.9 


28 


63 ±0 


04 


28 


33 ±0 


05 


28 


51 ±0 


08 






283.9 


29 


05 ±0 


05 


28 


61 ±0 


06 


29 


04itO 


13 






296.3 


29 


19 ±0 


06 


28 


78 ±0 


07 


30 


47 ±0 


42 






309.0 


29 


70 ± 


10 


29 


38 ±0 


12 












322.1 


29 


85 ±0 


11 


29 


04 ±0 


09 












335.6 


30 


44 ±0 


17 


29 


33 ±0 


11 












349.5 








31 


36 ±0 


55 
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Table 10 

SurfELce brightness profiles of NGC7716 



Radius 








Mi' 




[arcsec] 


[mag arcsec j 


[mag arcsec j 


r — 2i 

[mag arcsec J 


[mag arcsec J 


[mag arcsec j 


2.0 


on KQ -1- n m 


J-O.o^ _i_ ^J^\JlJ 


1 Q ni -\-f)f)A 


1 7 fin _|_ n 04 


1 7 Q4 -L n 04 

X 1 .tj^ _1_ \J.\J'-t 


5.9 


22.11 ± 0.01 


20.32 ± 0.01 


19.60 ±0.01 


19.17 ± 0.01 


18.97 ±0.01 


10.0 


22.63 ± 0.00 


21.00 ± 0.00 


20.32 ± 0.00 


19.90 ± 0.00 


19.73 ± 0.00 


14.1 


22.81 ± 0.01 


21.33 ± 0.01 


20.71 ± 0.01 


20.32 ± 0.00 


20.16 ±0.00 


18.4 


23.33 ± 0.01 


21.82 ± 0.01 


21.18 ±0.00 


20.78 ± 0.00 


20.64 ± 0.00 


22.8 


23.31 ±0.01 


21.96 ± 0.00 


21.42 ±0.00 


21.06 ± 0.00 


20.94 ± 0.00 


27.3 


23.59 ±0.01 


22.26 ± 0.00 


21.74 ±0.00 


21.40 ±0.00 


21.30 ±0.00 


31.9 


23.87 ± 0.01 


22.57 ±0.00 


22.09 ± 0.00 


21.75 ± 0.00 


21.67 ± 0.00 


36.7 


24.15 ± 0.01 


22.86 ± 0.00 


22.41 ± 0.00 


22.08 ± 0.00 


21.99 ±0.00 


41.7 


24.41 ±0.01 


23.19 ± 0.00 


22.77 ± 0.00 


22.45 ± 0.00 


22.38 ±0.00 


46.8 


24.87 ± 0.01 


23.58 ± 0.00 


23.16 ± 0.00 


22.86 ± 0.00 


22.80 ±0.01 


52.0 


25.26 ± 0.01 


23.99 ± 0.00 


23.56 ± 0.00 


23.25 ± 0.00 


23.20 ± 0.01 


57.4 


25.86 ± 0.02 


24.48 ± 0.00 


24.02 ± 0.00 


23.70 ± 0.00 


23.64 ± 0.01 


63.0 


26.25 ± 0.02 


24.88 ± 0.00 


24.39 ± 0.00 


24.07 ± 0.00 


24.00 ±0.01 


68.7 


26.66 ± 0.04 


25.30 ± 0.01 


24.80 ± 0.01 


24.49 ± 0.01 


24.42 ± 0.02 


74.6 


26.97 ± 0.05 


25.67 ±0.01 


25.16 ± 0.01 


24.83 ± 0.01 


24.74 ± 0.03 


80.7 


27.38 ± 0.06 


25.93 ± 0.01 


25.38 ± 0.01 


25.07 ± 0.01 


25.04 ±0.03 


87.0 


27.55 ± 0.08 


26.12 ± 0.01 


25.55 ± 0.01 


25.28 ± 0.01 


25.20 ± 0.04 


93.4 


27.73 ± 0.09 


26.31 ± 0.01 


25.77 ±0.01 


25.52 ±0.01 


25.40 ± 0.05 


100.1 


27.84 ± 0.09 


26.59 ± 0.01 


26.08 ± 0.01 


25.85 ±0.02 


25.73 ± 0.06 


106.9 


28.09 ± 0.10 


26.85 ± 0.01 


26.29 ± 0.01 


26.06 ± 0.02 


25.77 ± 0.05 


113.9 


28.36 ± 0.12 


27.03 ± 0.02 


26.47 ± 0.01 


26.26 ± 0.02 


26.02 ± 0.06 


121.2 


28.33 ±0.11 


27.01 ± 0.01 


26.47 ±0.01 


26.27 ±0.02 


26.16 ±0.07 


128.7 


28.20 ± 0.10 


27.01 ± 0.01 


26.44 ± 0.01 


26.33 ± 0.02 


26.06 ± 0.06 


136.4 


28.41 ± 0.11 


27.40 ± 0.02 


26.81 ± 0.02 


26.58 ± 0.02 


26.47 ± 0.08 


144.3 


29.24 ± 0.20 


27.94 ± 0.03 


27.34 ± 0.02 


27.20 ± 0.03 


26.95 ± 0.11 


152.5 


29.11 ± 0.18 


28.53 ± 0.04 


27.63 ± 0.03 


27.61 ± 0.04 


26.82 ±0.09 


160.9 


29.27 ± 0.20 


28.60 ± 0.05 


27.76 ± 0.03 


27.86 ± 0.06 


26.95 ±0.10 


169.6 


29.49 ± 0.22 


28.57 ±0.04 


27.86 ± 0.03 


27.72 ± 0.05 


27.05 ± 0.11 


178.5 




28.85 ± 0.05 


27.91 ± 0.03 


27.86 ± 0.05 


27.07 ±0.10 


187.7 




29.03 ± 0.06 


27.91 ± 0.03 


27.90 ± 0.05 


27.21 ±0.11 


197.2 




29.22 ± 0.07 


28.11 ±0.04 


27.87 ±0.05 


27.25 ±0.11 


206.9 




29.21 ± 0.06 


28.14 ±0.04 


27.95 ±0.05 


27.09 ±0.10 


217.0 




29.39 ± 0.07 


28.29 ± 0.04 


27.84 ± 0.04 


27.23 ± 0.10 


227.3 




29.47 ±0.08 


28.46 ± 0.04 


28.06 ± 0.05 


27.46 ±0.13 


238.0 






28 66 + n 06 




27 cjfi + n 14 


249.0 




30.41 ±0.18 


28.75 ±0.06 


28.46 ± 0.07 


27.42 ±0.12 


260.3 




29.68 ±0.09 


28.79 ±0.06 


28.46 ± 0.07 


27.77 ±0.16 


271.9 




29.69 ±0.09 


28.90 ±0.06 


28.55 ± 0.07 


27.57 ±0.13 


283.9 




30.13 ±0.12 


29.16 ±0.08 


28.98 ±0.10 


27.88 ±0.17 


296.3 




32.10 ±0.56 


30.14 ±0.16 


30.59 ± 0.36 


27.98 ±0.17 


309.0 






29.99 ±0.14 




28.93 ± 0.36 


322.1 






29.90 ±0.12 






335.6 






30.17 ±0.15 






349.5 






31.03 ±0.30 
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Table 11 

Surface brightness profiles of UGC02081 



Radius iXgi n^' Mi' Mz' 

[arcsec] [mag arcsec"^] [mag axcsec"^] [mag arcsec"^] [mag arcsec"^] 



2.0 


21 


71 


± 0.01 


21 


27 ± 0.01 


20 


89 ± 


02 


21 


02 


-I- 





02 


5.9 


22 


21 


±0.01 


21 


79 ± 0.01 


21 


42 ± 


01 


21 


58 


± 





01 


10.0 


22 


56 


± 0.00 


22 


18 ± 0.00 


21 


81 ± 


00 


21 


99 


± 





00 


14.1 


22 


73 


± 0.00 


22 


37 ± 0.00 


22 


02 ± 


00 


22 


19 


± 





00 


18.4 


22 


90 


±0.00 


22 


57 ± 0.00 


22 


22 ±0 


00 


22 


42 


± 





00 


22.8 


23 


11 


±0.00 


22 


78 ± 0.00 


22 


44 ±0 


00 


22 


65 


± 





00 


27.3 


23 


37 


± 0.00 


23 


04 ± 0.00 


22 


69 ±0 


00 


22 


90 


± 





00 


31.9 


23 


71 


± 0.00 


23 


37 ± 0.00 


23 


02 ± 


00 


23 


22 


± 





01 


36.7 


23 


93 


± 0.00 


23 


61 ± 0.00 


23 


27 ±0 


00 


23 


47 


± 





01 


41.7 


24 


17 


±0.00 


23 


85 ± 0.00 


23 


51 ±0 


00 


23 


73 


± 





01 


46.8 


24 


43 


± 0.00 


24 


12 ± 0.00 


23 


77 ±0 


00 


23 


98 


± 





01 


52.0 


24 


77 


± 0.00 


24 


45 ± 0.00 


24 


10 ±0 


00 


24 


30 


± 





01 


57.4 


25 


15 


± 0.00 


24 


81 ± 0.00 


24 


43 ± 


00 


24 


64 


± 





01 


63.0 


25 


45 


±0.00 


25 


14 ± 0.00 


24 


76 ±0 


01 


24 


99 


± 





02 


68.7 


25 


88 


± 0.01 


25 


56 ± 0.01 


25 


19 ±0 


01 


25 


44 


± 





02 


74.6 


26 


44 


± 0.01 


26 


09 ± 0.01 


25 


63 ± 


01 


25 


76 


± 





03 


80.7 


27 


12 


± 0.02 


26 


76 ± 0.02 


26 


20 ±0 


02 


26 


48 


± 





06 


87.0 


27 


81 


±0.03 


27 


47 ± 0.03 


26 


70 ±0 


02 


27 


24 


± 





11 


93.4 


28 


45 


±0.05 


28 


10 ± 0.05 


27 


13 ±0 


04 


27 


66 


± 





16 


100.1 


28 


58 


±0.05 


28 


49 ± 0.07 


27 


26 ±0 


04 


27 


51 


± 





14 


106.9 


29 


34 


± 0.10 


28 


80 ± 0.09 


27 


49 ± 


05 


27 


67 


± 





16 


113.9 


29 


78 


± 0.14 


29 


19 ± 0.12 


27 


92 ± 


07 


27 


71 







16 


121.2 


29 


94 


±0.16 


29 


17 ±0.11 


27 


65 ±0 


05 


27 


85 


± 





18 


128.7 


30 


59 


±0.25 


29 


26 ±0.10 


27 


77 ±0 


05 


28 


21 


± 





21 


136.4 








29 


37 ±0.12 


27 


92 ±0 


06 












144.3 








29 


82 ±0.16 


28 


05 ±0 


06 












152.5 








29 


81 ±0.15 


28 


20 ±0 


07 












160.9 








30 


60 ± 0.29 


28 


24 ±0 


07 












169.6 












28 


47 ±0 


08 












178.5 












28 


66 ±0 


10 












187.7 












29 


00 ±0 


13 












197.2 












29 


15 ±0 


14 












206.9 












29 


67 ±0 


21 
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Table 12 

SurfELce brightness profiles of UGC02311 



Radius 














[arcsec] 


[mag arcsec 


[mag arcsec"^] 


[mag arcsec"^] 


[mag arcsec ^ 


2.0 


19.95 di 0.03 


19.23 lb 0.03 


18.83 lb 0.03 


18 


57 ±0 


02 


o.y 


on Qc _i_ n ni 


on 91 -1- n ni 
zu.zi zn U.UI 


1 Q sn -1- n ni 
ly.ou m U.UI 


19 


57 ±0 


01 


inn 


91 c;q -I- n ni 
zi.oy m u.ui 


9n Q9 J_ n ni 

zu.yz zc U.UI 


90 -U n ni 
zv.oo m U.UI 


20 


34 ± 


01 


1/11 


oi QT -I- n ni 

ZL.o ( ± U.UI 


91 '39 _|_ n m 
Zl.oZ It U.UI 


9n -I- n Ai 
zu.yo it U.UI 


20 


81 ±0 


01 


10.41: 


OO yi? _|_ n ni 
ZZ.^i ± U.UI 


91 Q1 -U n m 
zi.yi it U.UI 


91 c^? -1- n ni 

Zl.O ( it U.UI 


21 


42 ±0 


01 


ZZ.o 


OQ 1 -1- n m 

Zo.LD ± U.UI 


zz.oo ± U.UI 


99 1 Q _i_ n m 

ZZ.Lo It U.UI 


21 


99 ±0 


01 




0*^ p;i -U n ni 

Zo.Ol- ± U.UI 


99 Qc _i_ n m 

ZZ. oO it U.UI 


99 Af\ J- n m 

ZZ.40 It U.UI 


22 


31 ± 


01 


oi.y 


OQ Qo _i_ n nn 
Zo.oZ ± U.UU 


9Q 91 _i_ n nn 
Zo.Zl ± U.UU 


99 Q1 _|_ n nn 
ZZ.ol it U.UU 


22 


68 ± 


01 


OO. ( 


OA 1 n -U n nn 

Z4.1U it U.UU 


9*3 _|_ n nn 
ZtD.O ( it U.UU 


9Q 91 _|_ n nn 
Zo.Zl it U.UU 


23 


11±0 


01 


41. / 


0/1 ^'i _i_ n ni 

Z4.00 it U.UI 


OA n/t -U n m 

Z^.XjH: It U.UI 


9Q '7n _i_ n m 

Zo. / U it U.UI 


23 


63 ±0 


01 


40.0 


oc; in-i-nni 

ZO.IU It U.UI 


9/1 f:l _|_ n m 
Z4.01 It U.UI 


9/1 9/1 _|_ n m 
Z4.Z4 It U.UI 


24 


18 ± 


01 


c^o n 


zo. / o it U.UI 


9c: 9Q _|_ n m 
zo.zo it U.UI 


9/1 Qfi _i_ n m 
z4.oD it U.UI 


24 


83 ±0 


02 


57.4 


26.59 ± O.Oi 


io.Uo ± U.UI 


25.63 ± 0.01 


25 


66 ±0 


04 


63.0 


27.32 lb 0.03 


26.74 lb 0.02 


26.42 ± 0.03 


26 


52 ±0 


09 


68.7 


27.82 lb 0.04 


27.32 lb 0.03 


26.81 ±0.04 


27 


16 ±0 


14 


74.6 


28.32 lb 0.05 


27.63 lb 0.04 


27.50 ±0.06 


27 


77 ±0 


20 


80.7 


28.67 lb 0.07 


27.90 lb 0.04 


27.71 ±0.07 








87.0 


28.95 lb 0.08 


28.19 lb 0.05 


28.02 ± 0.08 








93.4 


28.62 lb 0.06 


28.08 lb 0.05 


27.87 ± 0.07 








100.1 


29.62 lb 0.14 


28.76 lb 0.08 


29.25 ± 0.22 








106.9 


29.37 lb 0.11 


28.87 lb 0.09 










113.9 


29.89 lb 0.16 


29.65 lb 0.16 











